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FINAL  TECHNICAL  REPORT 

"COHERENT  COOPERATIVE  RADIATION  PROCESSES  IN  SOLIDS" 
(November  1,  1993  -  March  31,  1995) 


PROJECT  DESCRIPTION 

This  program  is  intended  to  explore  cooperative  interactions  among  rare  earth  dopants  in 
dielectric  solids,  with  upconversion  lasers  and  nonlinear  optical  applications  in  mind.  The 
focus  is  on  fundamental  studies  of  multi-atom  upconversion  and  avalanche  processes  rather 
than  multi-photon  upconversion  for  new  short  wavelength  lasers.  Because  of  significant 
potential  applications  in  the  areas  of  short  wavelength  sources,  optical  memories, 
communications,  and  displays  it  is  important  to  improve  our  understanding  of  cooperative 
dynamics  in  regard  to  their  optical  and  spatial  coherence,  temporal  evolution  and  relaxation 
mechanisms. 


FINAL  YEAR  SUMMARY 

In  the  final  year  of  this  three-year  project,  the  first  observations  of  ultrafast  dephasing  on 
cooperative  transitions  were  made  and  optical  switching  mediated  by  both  cooperative 
upconversion  and  avalanche  interactions  was  reported  for  the  first  time.  The  latter  work  in 
particular  represented  a  breakthrough  in  experimental  control  over  nonlinear  interactions  in 
dense  media,  revealing  that  excited  state  dipole-dipole  interactions  can  enhance  the  Lorentz 
local  field  in  dense  media  in  a  very  fundamental  way,  and  even  cause  intrinsic  bistability. 
The  coherent  transient  experiments  revealed  previously  unknown  dynamics  associated  with 
the  avalanche  effect,  important  for  a  full  understanding  of  the  behavior  of  avalanche 
upconversion  lasers. 

Near  the  beginning  of  the  third  year  of  this  project,  the  experimental  and  theoretical  effort 
turned  to  the  problem  of  careful  characterization  of  nonlinear  dynamics  in  cooperative 
interactions.  Two  main  tasks  were  undataken,  namely  direct  measurements  of  dephasing 
on  an  avalanche  transition  using  sub-picosecond  pulses,  and  the  observation  of  optical 
switching  and  instabilities  mediated  by  electromagnetic  coupling  between  atoms  in  dense 
media. 

On  the  topic  of  dephasing  it  was  shown  that  polarization  of  avalanche  transitions  in 
Tm:LiYF4  and  Tm:YAlC^  decayed  on  an  unexpectedly  fast  timescale.  In  YAIO3,  coherence 
decayed  in  60  ps  or  less  at  low  temperatures,  and  in  LiYF4  dephasing  was  too  rapid  to  be 
resolved  with  700  fs  pulses  below  temperatures  of  70K.  Detailed  measurements  of  T2 
versus  temperature  in  TmrYAlOs  were  analyzed  using  the  known  temperature  dependences 
of  various  dephasing  mechanisms  to  show  that  relaxation  proceeded  by  a  direct  phonon 
emission  process  induced  in  the  excited  state  by  the  avalanche  itself.  These  dynamics  have 
not  been  reported  previously,  and  strongly  affect  the  depletion  of  the  lower  avalanche  level. 
For  this  reason  our  findings  are  important  in  the  context  of  upconversion  lasers.  Similar 
considerations  should  apply  to  the  dynamics  responsible  for  fast  decay  in  Tm;LiYF4  but 
shorter  pulses  will  be  requued  in  future  measurements  to  verify  this. 

On  the  topic  of  optical  switching,  the  first  experiments  attempting  to  demonstrate  all-optical 
switching  mediated  by  avalanche  or  other  cooperative  interactions  were  performed.  A  self¬ 
pulsing  instability  was  generated  in  a  thin  disk  of  Tm:LiYF4  when  excited  at  the  avalanche 
transition  in  the  red  spectral  region  inside  an  enhancement  cavity.  Losses  in  the  external 
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cavity  were  minimized  by  working  with  a  100  pm  thick  sample  mounted  inside  at 
Brewster's  angle.  A  Drever-Pound  frequency-offset  locking  scheme  was  implemented  to 
permit  observations  at  arbitrary  detunings  of  the  cavity  from  laser  resonance,  and 
transmission  measurements  versus  time  were  then  made.  The  periodic  transmission 
maxima  observed  in  these  experiments  were  explained  as  arising  from  an  avalanche 
nonlinearity  opposed  by  a  competing  thermal  nonlinearity.  These  results  constitute  the  first 
reports  of  cavity -enhanced  switching  in  an  avalanche  medium. 

In  parallel  with  these  developments,  it  was  also  discovered  that  intrinsic  bistability  resulted 
from  cooperative  emission  in  crystals  of  Yb3+:CsY2Br9  at  low  temperatures.  Green 
upconversion  luminescence  intensity  was  measured  versus  temperature  and  versus  input 
power.  In  this  way  abrupt  steps  in  emission  were  noticed  which  were  qualitatively 
reproduced  by  rate  equations  provided  specific  assumptions  were  made  about  temperature 
dependences  of  upconversion  and  migration  coefficients.  However  this  phenomenological 
approach  gave  no  suggestions  regarding  the  basic  origin  of  the  observed  switching  and 
considerable  effort  was  expended  developing  a  complete  density  matrix  theory  of  this 
effect.  One  of  the  greatest  successes  of  this  program  to  date  was  the  demonstration  that  a 
density  matrix  treatment,  and  more  specificdly  a  correct  treatment  of  dephasing 
interactions,  was  necessary  to  obtain  complete  agreement  with  experiments.  'ITie  density 
matrix  theory  revealed  that  fundamental  roles  of  upconversion  and  migration  interactions  in 
the  excited  state  were  opposite.  Resonant  migration  tended  to  suppress  optical  hysteresis, 
whereas  upconversion  generated  it. 

The  data  also  confirmed  that  the  mechanism  responsible  for  intensity-dependent  optical 
switching  was  an  enhancement  by  upconversion  of  the  implicit  nonlinearity  in  optical 
response  of  all  media  caused  by  the  difference  between  the  local  field  and  the  input  electric 
field  of  the  light.  That  is,  our  experiments  established  for  the  first  time  that  excited  state 
dipole-dipole  interactions  can  have  a  profound  effect  on  the  Lorentz  field  of  all  media, 
particularly  dense  media. 
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PUBLICATIONS  TO  DATE 

1 .  H.  Ni  and  S.  C.  Rand,  "Avalanche  Upconversion  in  Tm: YALO,"  Opt.  Lett.  16, 1424 
(1991). 

2.  D.  Redman,  S.  Brown,  R.  Sands  and  S.  C.  Rand,  "Spin  Dynamics  and  Electronic 
States  of  N-V  Centers  in  Diamond,"  Phys.  Rev.  Lett.  62,  3420  (1991). 

3.  D.  Redman,  S.  Brown  and  S.  C.  Rand,  "Origin  of  Persistent  Hole-Burning  of  N-V 
Centers  in  Diamond,"  J.  Opt.  Soc.  Am.  B  2,  No.  5,  (1992). 

4.  D.  Redman,  Q.  Shu,  S.  W.  Brown,  A.  Lenef,  Y.  Liu,  J.  Whitaker,  S.  C.  Rand,  S. 
Satoh,  K.  Tsuji  and  S.  Yazu,  "Electronic  Structure  of  N-V  Centers  and  Terahertz 
Spectroscopy  of  Diamond,"  AIP  Proceedings  Volume  (1992). 

5.  D.  Redman,  Q.  Shu,  A.  Lenef,  and  S.  C.  Rand,  "Two-Beam  Coupling  by  Nitrogen- 
Vacancy  Centers  in  Diamond,"  Opt.  Lett.  12, 175  (1992). 

6.  P.  Xie  and  S.  C.  Rand,  "Astigmatically-Compensated,  High  Gain  Cooperative 
Upconversion  Laser,"  Appl.  Phys.  Lett.  60,  3084(1992). 

7.  P.  Xie  and  S.  C.  Rand,  "Visible  Cooperative  Upconversion  Laser,"  Opt.  Lett.  17, 
(1992). 

8.  H.  Hi,  Q.  Shu,  and  S.  C.  Rand,  "Observation  of  the  Cooperative  Kerr  Effect,"  Phys. 
Rev.  Lett,  (in  preparation). 

9.  H.  Ni  and  S.C.  Rand,  "Avalanche  Phase  Conjugation",  Opt.  Lett.  12,  No.  17(1992). 

10.  P.  Xie  and  S.C.  Rand,  "Cw  Mode-locked  Visible  Upconversion  Laser",  Opt.  Lett. 
12,  No.  16(1992). 

11.  P.  Xie  and  S.C.  Rand,  "Nonlinear  Dynamics  of  Cooperative  Upconversion", 

J.O.S.A.  B  11.  901(1994). 

12.  P.  Xie  and  S.C.  Rand,  "Continuous-wave,  four-fold  Upconversion  Laser",  Appl. 

Phys.  Lett.  62,  3125(1993). 

13.  Q.  Shu,  H.  Ni  and  S.C.  Rand,  "Avalanche  Two-beam  Coupling",  Optics  Letters 

(submitted). 

14.  Q.  Shu,  S.W.  Brown,  J.  Bretz,  S.C.  Rand,  L.-Q.  Zhang  and  M.  Bretz,  "Intersystem 
Crossing  of  Ceo-  Ullt^igh  Resolution  Differential  Transmission  Spectroscopy", 

Opt.  Lett,  (submitted). 

15.  A.  Lenef,  S.W.  Brown,  and  S.C.  Rand,  "Implications  of  Zeeman  Coherence  and 
Quantum  Beats  in  Ultrafast  Photon  Echoes  of  N-V  Centers  in  Diamond",  Phys.  Rev. 
Lett,  (in  preparation). 


16.  A.  Lenef,  D.  Kreysar,  K.  Obermyer  and  S.C.  Rand,  "Laser-induced  Collisional 
Avalanche  in  Atomic  Cesium",  Phys.  Rev.  Lett,  (submitted). 

17.  S.C.  Rand,  A.  Lenef,  and  S.W.  Brown,  "Zeeman  coherence  and  quantum  beats  in 
ultrafast  photon  echoes  of  N-V  centers  in  diamond",  J.  Luminescence  60&61. 
739(1994). 

18.  A.  Lenef  and  S.C.  Rand,  "Tunneling  Dynamics  of  Squeezed  States  in  a  Potential 
Well",  Phys.  Rev.  A  49.  32(1994). 

19.  M.  Hehlen,  H.U.  Gudel,  J.  Rai,  S.  Rai  Q.  Shu  and  S.  Rand,  "Cooperative 
Bistability  in  Dense,  Excited  Atomic  Systems",  Phys.  Rev.  Lett.  22, 1 103(1994). 

20.  M.  Hehlen,  H.  Gudel,  Q.  Shu  and  S.C.  Rand,  Cooperative  Bistable  Dynamics  in 
Yb3+:CsY2Br9",  Phys.  Rev.  B  (in  preparation). 

21 .  H.  Ni  and  S.C.  Rand,  "Avalanche  Beam  Propagation",  Opt.  Lett,  (in  preparation). 

22.  S.C.  Rand  (invited).  Diamond  Lasers",  in  Properties  and  Growth  of  Diamond,  ed. 

G.  Davies,  linstitution  of  Electrical  Engineers  (Br.),  London,  United  Kingdom, 

1994,  pp.  235-239. 

23.  I.  McMichael,  R.  Saxena,  T.Y.  Chang,  Q.  Shu,  S.C.  Rand,  J.  Chen  and  H.  Tuller, 
High  Gain  Nondegenerate  Two-wave  Mixing  in  CriYAlOa",  Opt.  Lett,  (submitted). 

GRADUATE  STUDENTS  SUPPORTED: 

(Fall  '91  -  Fair93) 

TOTAL  NO.  GRAD  STUDENTS  WHO  RECEIVED  AFOSR  SUPPORT=  10 

H.  Ni-Physics-Ph.D.  Degree  1994 

A.  Lenef- Applied  Physics-Ph.D.  Degree  1993 
S.  Brown-Applied  Physics-Ph.D.  Degree  1994 
Qize  Shu-Applied  Physics-Candidate 
Ping  Xie- Applied  Physics-Ph.D.  Degree  1992 
D.  Kreysar-Applied  Physics-M.Sc.  Degree  1993 
S.  Kidner-Physics-Ph.  D.  1994  (co-chair) 

Lih-Mei  Yang-Physics-Candidate 
Dave  Redman  -EECS-Ph.D.  Degree  1991 
K.  Obermyer-Physics-Pre-candidate 

THESES: 

I.  D.A.  Redman,  "Electronic  Structure  of  the  Nitrogen- Vacancy  Color  Center  in 
Diamond",  Ph.  D.  Dissertation,  University  of  Michigan,  1991. 

2.  P.  Xie,  "Continuous- wave  Cooperative  Upconversion  Lasers",  Ph.  D.  Dissertation, 
University  of  Michigan,  1992. 

3.  A.  Lenef,  "Ultrafast  Photon  Echo  Spectroscopy  of  Luminescent  Centers:  N-V:Diamond 
and  Tm^+:YA103",  Ph.  D.  Dissertation,  University  of  Michigan,  1993. 

4.  H.  Ni,  "Avalanche  Upconversion  in  Tm:LiYF4,  Tm:YA103  and  Tm:YAG",  Ph.D. 
Dissertation,  University  of  Michigan,  1994. 


5.  S.W.  Brown,  "Electronic  Structure  of  Luminescent  Centers  in  Wide  Bandgap 
Semiconductors",  Ph.  D.  Dissertation,  University  of  Michigan,  1994. 

6.  S.  Kidner,  "Ion-assisted,  Sputter  Deposition  and  Structural  Characterization  of  Cubic 
Boron  Nitride",  Ph.  D.  Dissertation,  University  of  Michigan,  1994. 

PATENTS: 

Continuous- wave,  pair-pumped  laser  (Patent  No.5, 116,437). 

Continuous-wave,  trio-pumped  laser  (Patent  No.  5,038,358). 

Four-fold  and  Higher  C)rder,  Continuous- wave  Upconversion  Lasers  (Patent  No.5,289,481). 
Mode-locked  Upconversion  Laser  Source  (Patent  No.5,27 1,058). 

INVITED  CONFERENCE  PRESENTATIONS: 

S.C.  Rand(invited),  "Cooperative  Upconversion  Lasers  and  Nonlinear  Dynamics",  OSA 
Annual  Meeting,  Albuquerque,  New  Mexico,  Sept.  20-25, 1992,  paper  M12. 

S.C.  Rand  (invited),  "A  Hydrogenic  Model  of  the  N-V  Center  in  Diamond",  Annual 
Meeting  of  the  American  Physical  Society,  Indianapolis,  March  16-20, 1992,  paper  029-1. 

S.C.  Rand  (invited),  "Upconversion  Lasers",  French-US  Symposium  on  Laser-Matter 
Interactions,  University  of  Michigan,  Ann  Arbor,  Michigan,  April  29-30, 1993. 

S.C.  Rand  (invited),  "Spin  Dynamics,Quantum  Beats  and  Electronic  States  of  Luminescent 
Centers  in  Diamond",  Ann.  Mtg.  of  Fed.  Analyt.  Chem.  and  Spectr.  Soc.  (FACSS’93), 
Detroit,  Michigan,  Oct.  17-22, 1993. 

OTHER  CONFERENCE  PRESENTATIONS: 

P.  Xie  and  S.C.  Rand,  "Continuous-wave,  pair-pumped  Laser",  in  Advanced  Solid  State 
Lasers,  Vol.  6,  eds.  H.  Jenssen  and  G.  Dube,  Optic^  Society  of  America,  1991,  pp.  190- 
191. 

D.A.  Redman,  Q.  Shu,  S.W.  Brown,  A.  Lenef,  Y.  Liu,  J.  Whitaker,  S.  Rand,  S.  Staoh, 

K.  Tsuji  and  S.  Yazu,  "Electronic  States  of  N-V  Centers  and  Terahertz  Spectroscopy  of 
Diamond",  Proc.  of  Mat.  Res.  Soc.,  Fall  Meeting,  Boston,  Massachusetts,  November 
1991,  paper  G1.5. 

P.  Xie  and  S.C.  Rand,  "Cooperative  Nonlinear  Dynamics  of  the  2.8  |j.m  pair-pumped 
erbium  laser,  Conference  on  Lasers  and  Electro-optics  (CLEO'92),  Anaheim,  California, 
May  10-15, 1992,  paper  CFE4. 

Q.  Shu,  S.  Brown,  A.  Lenef,  D.  Redman,  and  S.C.  Rand,  "Nonlinear  Dynamics  and 
Electronic  States  of  N-V  Centers  in  Diamond",  Ann.  Mtg.  American  Physical  Society, 
Indianapolis,  Indiana,  March  16-20, 1992,  paper  126  8. 

S.C.  Rand,  D.  Redman,  and  S.W.  Brown,  "Origin  of  Persistent  Hole-burning  of  N-V 
Centers  in  Diamond",  American  Physical  Society,  March  Meeting,  Indianapolis,  Indiana, 
March  16-20, 1992. 


6 


A.  Lenef  and  S.C.  Rand,  "Quantum  Beats  and  Spin  Coherences  in  N-V  Centers  in 
Diamond",  Ann.  Mtg.  American  Physical  Society,  Indianapolis,  Indiana,  March  16-20, 

1992,  paper  126  9. 

S.C.  Rand,  D.A.  Redman  and  S.W.  Brown,  "Origin  of  Persistent  Hole-burning  of  N-V 
Centers  in  Diamond",  Ann.  Mtg.  American  Physical  Society,  Indianapolis,  Indiana,  March 
16-20,  1992,  paper  126  10. 

H.  Ni  and  S.C.  Rand,  "Spatial  Ring  Formation  in  Avalanche  Optical  Nonlinearities",  Ann. 
Mtg.  American  Physical  Society,  Indianapolis,  Indiana,  March  16-20, 1992,  paper  018  - 
11. 

S.C.  Rand,  "Nonlinear  Dynamics  of  Upconversion  Lasers",  Program  in  Nonlinear 
Studies,  Poster  Conference,  University  of  Michigan,  Ann  Arbor,  Michigan,  January  19, 

1993. 

I.  McMichael,  R.  Saxena,  T.Y.  Chang,  R.  Neurgaonkar,  and  S.C.  Rand,  "High  gain, 
non-degenerate  Two-Wave  Mixing  in  Cr:YA103",  Conf.  on  Lasers  and  Electro-optics 
(CLEO'93),  Baltimore,  Maryland,  May  2-7,  1993,  paper  CThS34. 

H.  Ni,  Q.  Shu  and  S.C.  Rand,  "Avalanche  Phase  Conjugation",  Qu.  Elect,  and  Laser 
Science  Conf.  (QELS'93),  Baltimore,  Maryland,  May  2-7, 1993,  paper  QThD5. 

A.  Lenef,  S.  Brown,  and  S.C.  Rand,  "Quantum  Beats  in  Diamond",  Qu.  Elect,  and  Laser 
Science  Conf.  (QELS'93),  Baltimore,  Maryland,  May  2-7, 1993,  paper  QThH35. 

S.C.  Rand,  A.  Lenef,  and  S.W.  Brown,  "Zeeman  Coherence  and  Quantum  Beats  in 
Ultrafast  Photon  Echoes  of  N-V  Centers  in  Diamond",  Int.  Conf.  on  Lumin.  (ICL'93), 
Storrs,  Connecticut,  Aug.  9-13,  1993. 

S.C.  Rand,  A.  Lenef,  and  S.W.  Brown,  "Excited  State  Structure  and  Dephasing  of  Point 
Defects  in  Widegap  Semiconductors:  Ultrafast  Four-Wave  Mixing  Spectroscopy  of  N-V 
Centers  in  Diamond",  Tech.  Conf.  2041  of  Soc.  Photo-Optical  Ristrum.  Eng.  (S.P.I.E.), 
Quebec  City,  Canada,  Aug.  16-20,  1993,  paper  2041-18. 

S.C.  Rand,  H.  Ni,  Q.  Shu  and  P.  Xie,  "Mode-locking  Applications  of  Nonlinear 
Cooperative  Interactions",  Tech.  Conf.  2041  of  Soc.  Photo-Optical  Instrum.  Eng. 
(S.P.I.E.),  Quebec  City,  Canada,  Aug.  16-20,  1993,  paper  2041-18. 

A.  Lenef,  D.  Kreysar,  and  S.C.  Rand,  "Collisional  Avalanche  in  Atomic  Cesium",  Opt. 
Soc.  of  America  (OSA’93),  Toronto,  Canada,  Oct.  3-8, 1993,  paper  MTT7. 

A.  Lenef  and  S.C.  Rand,  "Tunneling  Dynamics  of  Squeezed  States  in  a  Potential  Well", 
Opt.  Soc.  of  America  (OSA’93),  Toronto,  Canada,  Oct.  3-8, 1993,  paper  ThL3. 

S.C.  Rand  and  P.  Xie,  "Intrinsic  Instabilities  of  Cooperative  Upconversion  Dynamics  in 
Lasers",  Opt.  Soc.  of  America  (OSA'93),  Toronto,  Canada,  Oct.  3-8, 1993,  paper  WSS5. 

S.C.  Rand,  A.  Lenef  and  D.  Kreysar,  "Atomic  Collisional  Avalanche",  Eleventh 
International  Conference  on  Laser  Spectroscopy,  Hot  Springs,  Virginia,  June  13-18, 1993, 
paper  MP3. 


M.  Hehlen,  H.U.  Gudel,  J.  Rai,  S.  Rai  and  S.C.  Rand,  "Cooperative  Bistability  in  Dense, 
Excited  Atomic  Sytems",  Int.  Qu.  Elect.  Conf.  (IQEC94),  Anaheim,  California,  May  8- 
13, 1994,  paper  QME8. 

A.  Lenef  and  S.C.  Rand,  "Photon  Echo  Measurements  of  Optical  Dephasing  on  a  Thulium 
Avalanche  Transition",  Int.  Qu.  Elect.  Conf.  (IQEC'94),  Anaheim,  California,  May  8-13, 
1994,  paper  QMC3. 

H.  Ni  and  S.C.  Rand,  "Temperature  Dependence  of  the  avalanche  threshold  and  the  cross 
relaxation  rate  in  Tm:LiYF4",  Int.  Qu.  Elect.  Conf  (IQEC'94),  Anaheim,  California,  May 
8-13,  1994,  paper  QTu05. 

I.  McMichael,  R.  Saxena,  T.Y.  Chang,  Q.  Shu,  S.C.  Rand,  J.  Chen,  and  H.  Tuller,  "high 
gain  Non-degenerate  Two-Wave  Mixing  in  Cr:  YAIO3",  OSA  Conf.  on  Nonlinear  Optics", 
Waikoloa,  Hawaii,  July  25-29, 1994,  paper  WAS. 


INVITED  COLLOQUIA: 

S.C.  Rand,  "Coherent  Cooperative  Radiation  Phenomena  in  Solids",  General  Colloquium, 
Dept,  of  Physics,  Alabama  A&M,  Huntsville,  Alabama,  June  16, 1992. 

S.C.  Rand,  "Ultrafast  and  Ultrahigh  Resolution  Spectroscopy  of  Diamond",  General 
Colloquium,  Dept,  of  Physics,  Alabama  A&M,  Huntsville,  Alabama,  June  15, 1992. 

S.C.  Rand,  "Cooperative  Nonlinear  Dynamics:  Upconversion,  Bistability,  and 
Delocalization  in  Rare  Earth  Systems",  General  Colloquium,  Dept,  of  Physics,  Wayne 
State  University,  November  1993. 

AWARDS: 

1993-  The  thesis  of  Ping  Xie  on  "Continuous- wave  Cooperative  Upconversion  Lasers" 
was  awarded  the  University  of  Michigan  Distinguished  Dissertation  Award. 

1993-  The  thesis  of  Ping  Xie  was  awarded  the  Terwilliger  Prize  by  the  Department  of 
Physics  of  the  University  of  Michigan. 

1994-  Rackham  Dissertation  Fellowship  for  Hui  Ni,  awarded  by  the  University  of 
Michigan. 

1994-PI  named  Fulbright  Scholar  by  the  Fulbright  Foundation 

1994-PI  named  a  Fellow  of  the  Joint  Institute  of  Laboratory  Astrophysics  (University  of 
Colorado  and  JILA). 


SECTION  2 


REPRINTS 
(1993-1994  only) 


a  reprint  from  Journal  of  the  Optical  Society  of  America  B 


901 


Nonlinear  dynamics  of  cooperative  upconversion 


P.  Xie*  and  S,  C.  Rand 

Division  of  Applied  Physics,  1049  Randall  Laboratory,  University  of  Wchigan,  Ann  Arbor,  Michigan  48109’!  120 


Received  July  15,  1993;  revised  manuscript  received  November  29,  1993 

We  describe  several  new  aspects  of  light -matter  interactions  for  solids  in  which  interatomic  coupling  of 
impurity  atoms  plays  a  dominant  role  in  population  dynamics.  We  explore  the  implications  of  spatial 
coherence  in  such  multiatom  interactions  by  introducing  a  density-matrix  theory  of  cooperative  upconversion, 
focusing  on  pair  systems  for  which  analytic  results  can  be  obtained.  We  predict  population  pulsations  in 
coherent  cooperative  upconversion,  enhanced  quantum  efficiency,  enhanced  energy  transfer,  and  pair-mediated 
instabilities,  not  only  in  cooperative  upconversion  media  without  external  cavities  but  in  upconversion  lasers 
and  conventional  lasers  in  highly  doped  solids  as  well.  These  predictions  are  compared  with  rate  equation 
solutions  and  observations  in  lasers  with  inversions  siistained  by  cooperative  processes,  particularly  the 
2.8- /xm  Er  laser.  Rate  equations  fail  to  predict  the  observed  steady-state  instabilities  of  this  laser,  which 
are  well  reproduced  by  density-matrix  theory,  furnishing  evidence  of  weak  coherent  delocalizations  in  a  rare- 
earth  system. 


1.  INTRODUCTION 

Upconversion  and  energy-transfer  processes  that  occur 
in  solids  as  a  result  of  nonlinear  population  dynamics, 
rather  than  parametric  processes,  have  been  studied  for 
many  decades. Two  classes  of  nonparametric  up¬ 
conversion  that  differ  in  the  mechanism  responsible  for 
converting  long-wavelength  radiation  incident  upon  a 
sample  into  short-wavelength  emission  are  generally  ac¬ 
knowledged.  In  the  first,  sequential  absorption  of  inci¬ 
dent  photons  by  progressively  higher  energy  levels  of  a 
single  atom  leads  to  emission  from  high-lying  upconverted 
states,  broadly  designated  as  upconversion  emission.  In 
the  second,  each  atom  absorbs  at  most  one  photon,  and 
energy  transfer  between  atoms  is  required  for  upconver¬ 
sion  emission.  The  latter  t)q>e  of  process  may  proceed  in 
one  of  three  ways.  In  a  cooperative  process  the  pooled 
energy  of  several  coupled  atoms  can  become  localized  on 
one  particular  atom,  which  then  emits  light  from  a  high- 
lying,  upconverted  state.^  Or  energy  may  be  transferred 
in  a  sequence  of  discrete  transfer  steps. ^  Another  pos¬ 
sibility  is  an  avalanche  process  in  which  upconversion 
results  from  field-induced  cross-relaxation  dynamics. 

In  this  paper  we  confine  our  discussion  to  the  cooperative 
upconversion  process,  highlighting  novel  experimental 
and  theoretical  d3mamics  that  can  result  from  spatial  co¬ 
herence  among  the  coupled  atoms  or  from  the  interaction 
of  coupled  atoms  with  radiation  fields  in  optical  cavities. 

Early  research  on  upconversion  was  motivated  by 
the  need  to  overcome  the  poor  sensitivity  of  infrared 
detectors.®  More  recent  research  on  upconversion  has 
been  stimulated  by  the  realization  that  upconversion 
lasers  have  potential  as  practical,  solid-state  sources  of 
visible  and  ultraviolet  light. However,  to  date,  little 
discussion  has  been  given  to  certain  fundamental  as¬ 
pects  of  the  upconversion  processes  themselves,  such  as 
basic  limits  to  achievable  degrees  of  upconversion  and  the 
inherent  stability  of  these  highly  nonlinear  interactions 
between  light  and  matter. 

In  this  paper  we  consider  these  topics,  primarily  in 
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the  context  of  new  experiments  with  the  pair-pumped 
Er  laser,  a  laser  with  an  inversion  sustained  purely  by 
pair  upconversion.^^  General  predictions  of  enhanced 
quantum  efficiency,  enhanced  energy  transfer,  popu¬ 
lation  pulsations,  and  instabilities  mediated  by  coop¬ 
erative  dynamics  are  compared  with  observations  of 
nonlinear  dynamics  in  twofold,  threefold,  and  fourfold 
upconversion  processes. 

2.  THEORY 

A.  Enhanced  Quantum  Efficiency 
The  general  scheme  of  cooperative  upconversion  is  illus¬ 
trated  in  Fig.  1.  Two  atoms  are  depicted  in  an  excited 
state  undergoing  a  pair  transition  in  which  one  atom  is 
promoted  to  an  upper  state  while  the  other  descends  to 
the  ground  state.  Multipole  or  exchange  coupling  is  pre¬ 
sumed  to  furnish  the  interaction  responsible  for  this  pair 
transition,  mediating  evolution  of  the  system  in  a  man¬ 
ner  somewhat  analogous  to  coupled  harmonic-oscillator 
dynamics.  The  energy  transferred  to  the  second  atom 
rarely  matches  the  precise  transition  energy  required  for 
gaining  access  to  a  real  final  state,  but  the  energy  defect 
can  be  accommodated  in  solids  by  the  emission  or  the  ab¬ 
sorption  of  phonons.  In  solids  doped  heavily  with  rare 
earth  or  transition-metal  impurities,  transitions  of  this 
kind  are  common.^ 

A  schematic  representation  of  excitation  and  decay 
pathways  involved  in  cooperative  upconversion  is  given 
in  Fig.  2.  Excitation  is  initiated  by  a  pumping  process 
at  wavelength  Am,  which  populates  level  |3)  with  effi¬ 
ciency  T/o.  The  system  subsequently  cascades  to  level 
|1>  in  steps  with  individual  branching  ratios  77,.  Decay 
to  the  groimd  state  |0)  may  then  occur,  or  a  cooperative 
transition  of  the  t3q>e  indicated  in  Fig.  1  may  take  place. 
The  manner  in  which  cooperative  upconversion  nt'fits 
the  overall  emission  at  Aout  is  not  immediatelv  .t  ^ious 
because  the  maximum  branching  ratio  for  a  pair 

upconversion  process  is  0.5.  This  seems  1ov\  ir  r.  the 
outset.  However,  it  can  readily  be  appreci.it* '^nm 
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Fig.  1.  Dynamics  of  cooperative  upconversion  in  a  system  con¬ 
sisting  of  two  four-level  atoms  initially  occupying  state  ll>,  as 
indicated  by  filled  circles.  The  system  reliro  by  promoting 
one  atom  to  level  |3)  as  the  second  decays  to  the  ground  state 
(curved  arrows).  In  this  fashion  upconversion  fluorescence  firom 
energy  levels  higher  than  the  initial  state  becomes  possible 
at  several  wavelengths  (strai^t  arrows).  The  carrespondence 
between  levels  of  the  model  and  Er  is  as  follows:  |0> 

ID  --  12)  -  13)  - 


Fig.  2.  Recycling  by  cooperative  upoonveraion.  The  ™tiel  ex¬ 
citation  of  level  |3),  either  by  resonant  optical  excitation  (aby  rn 
cooperative  process,  occurs  with  quantum  effideaaqr  This 
is  followed  by  a  cascade  of  nonraduitive  decays  and  emissioa 
processes  with  individual  branching  ratios  of  y^.  .  States  llaX 
and  lib)  are  different  Staric  components  of  le^  ID.  Finally 
a  long-lived  level  is  reached  in  which  coopesatwa  tqmmversioiir 
can  occur  spontaneously  witiwat  finther  abserptfam  d  photons 
fiom  the  pump  field.  Tbs  mnnant  exeitedkMta  population  can 
be  recycled  to  level  |3)  andjlwhtftkmal  jdiotoins  may  be  emittod* 
thereby  enhancing  quantiMd^fdBNAiKyf.  The  radiative  transi^m 
|2)  |lb)  is  the  focus  of  qgppfuHi-eflflckff^  omsideraticms. 

Fig.  2  that  excited-stata  population  can  bd  recyded 
through  the  emission  channel  if  spmtatneoua  upomver-^ 


Sion  occurs. 


To  quantify  this  concept  and  to  dioer  that  foantbrn 
efficiency  of  optical  exaiamm  can  be  snhaacedhy 
enoe  of  cooperative  dynamics,  we  asedme  that  tBe^too^pers^ 
tive  upconversion  pioeesa  can  repeat  itsebfitti^^ 
numb^  of  sites  at  which  upconverskm  occwa  is  tabfin 
to  be  small  compared  with  the  number  of  ezdted  atoms; 
This  ensures  that  the  reservoir  of  state  |1)  ions  surround¬ 
ing  upconversion  sites  remains  undepleted  throu^  many 
upconversion  cycles  because  of  rapid  energy  migration 
within  the  reservoir.  Ground-state  A  atoms  at  upconver- 
sion  sites  are  continually  reexdted  to  state  |1>  by  means 


of  resonant  energy  transfer  firom  the  reservoir,  restoring 
the  initial  condition  shown  in  Pig.  1  without  requiring  the 
absorption  of  additicmal  photons.  For  very  large  n  the 
effective  quantum  efficiency  then  becomes 

T/eff  =*  ^0’7i172[1  +  VbVmvP^^^ViV2 

+  VzVxxpP^^^mV2V3Vu^^^^VlV2  +  •••] 


X  V0VlV2iVBVxxpVlV2T 

n-0 


V0V1V2  . 

1  -  viV2Vinup 


Here  factors  P*"'  are  introduced  to  account  explicitly  for 
the  probability  of  A's  being  reexdted  by  energy  migration 
dnrilig  the  (^de.  Both  atoms  A  and  B  must  be  in 
state  |1>  for  upconversion  to  occur  with  branching  ratio 
^up*  Howevmr,  P^*  may  be  set  equal  to  unity  without  a 
significant  loss  of  ganerality  [P<^  •-  P<®  «...  P<*>  —  IJ 
For  a  fixed  initial  number  of  exdted  atoms,  P^”’  decreases 
as  the  reservdr  deplete^' it  is  an  eni^ent  approxi¬ 
mation  to  assume  that  P^is  independent  of  n  when  the 
density  of  upconveraios  attes  is  compared  with  the 
densi^  ci  exdted  ^om^  and  stelai^-state  iSumination 
is  assumed:  Wedm  omiqpm.th«  eflbctive  effideuiy  ii^r 
with  that  expected  without'piqEpdfdfim  reqrdii^hy  d^jfr 
ing  an  enhancement  fimtcr  Jf,  givan  1^  -  \ 


It  is  immediately  sqrparent  fidm  quanttuh 

effidendes  exceeding  nnify^  an  poasiIds.  '  TD  see  this, 
one  can  ignm  decay  fiwm  ^  upeonvaiBfam  stats  |la) 
by  setting  174  -  0  and  ^vahis  q.,  -•  (l/m)  -  i»4  “  IM 
into  E!q.  (2).  m  is  Qm'inadl^  oF  siiolns  partidpating  in 
the  upconveraion  pronm  ita^  Then,  for  example,  a 
quantum  ^Sdency  of’B  ia  expected^wfasn  branching  ra¬ 
tios  for  decay  and  emisitfah  in  equal  unify, 

and  is  taken  tn  b%VA  -  l^’  hfocanum  value  ap- 
propriato  for  pair  .upconvyafoBL^^  M  of  great  sig- 
nificanoe  because  pcpdtsHiai  also  yields  an 

enhanced  energy  convetskta  effidsdii^-  ||nn  by 

' '  r  O) 


FromEq.  (3)itisdsarfis^ftarn^4ddiiqps9r^veupcon- 
verskm  greatly  exceeding  the  nafuitf  decay  rate  of  state 
2,  energy  efficiency  eeaace  to  d^end  aonitlr^  on  the  ac¬ 
tual  vdtta  of  ffienpcoavwaion  wkyitiialf  Or-  <»  the  upcon- 
vermnbraaddi^ftdia'S^iiblia^toranhanoed  value. 

bdde  irreepee- 
axeinstrue- 


ty 

.anlltiiiil^itoFii 
ttUDfcfittOBODB; 

KiilS 

Ain/(m  -  1)  hiE 


1,  ftunished 
die  maximum 
^^Jbecauaa 
raadiea  ffie 
incycBng  tte 
ftum 
th*  ^ 
H^as  ahcM  as 
ar(w  -  1),  the 


energy  efficiency  can 

Second^  wlien  'a  .  aito^|ei^>i^v^  is  used  to 
pump  upper  levst  tj^bijpnwi^lpl^l^,  conventional  laser- 
pumping  sdieiMS'^ciBdliYitifitbljlji^^biljbfbwitnent  of  quantum 
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efficiency  can  still  occur,  provided  that  Ajn/Aout  ^  (m  - 
l)/m,  constraining  downward  and  upward  transitions  to 
conserve  energy  in  the  cooperative  process.  We  again 
find  that  the  maximum  conversion  efficiency  {r}E*)max  ^ 
{mAin/[(m  -  l)Aout]}nmx  =  1  excecds  the  limiting  value  set 
by  =  (Ain/Aout)mM  =  (m  “  l)/m  for  the  same  emis¬ 

sion  process  in  the  absence  of  recycling.  Energy  conver¬ 
sion  efficiency  t}e*  never  exceeds  unity  because  of  the 
restrictions  on  Ain/Aout  necessary  for  energy  transfer  to 
occur  between  excited  states. 

B.  Population  Pulsations 

Modifications  of  cooperative  dynamics  are  to  be  expected 
in  optical  cavities.  This  may  happen,  for  example,  if 
stimulated  emission  follows  cooperative  upconversion. 
For  a  stimulated  emission  process  that  repopulates  cou¬ 
pled  levels  while  depleting  the  gain,  a  sudden  increase 
in  the  cooperative  upconversion  rate  follows  gain  quench¬ 
ing.  The  abruptly  increased  rate  of  upconversion  tends 
to  restore  inversion,  and  a  relaxation  oscillation  can  take 
place. 

Rate  equation  calciilations  indicate  that  steady  re¬ 
laxation  oscillations  cannot  be  sustained  by  coopera¬ 
tive  d3mamics,  a  result  considered  in  more  detail  in 
Subsection  2.C.  However,  transient  oscillations  can  ap¬ 
pear  during  the  buildup  or  the  decay  of  light  in  the  cavity. 
We  may  show  this  by  solving  the  rate  equations  for  popu¬ 
lations  Ni  in  each  state  (i  =  0, 1, . . . ,  3)  together  with  the 
equation  for  cavity  photon  density  Q: 


In  these  expressions  yij  is  the  decay  constant  between 
levels  i  and  J,  yt  is  the  total  decay  rate  for  level  t,  Bij 
is  the  Einstein  coefficient  for  induced  absorption  on  the 
i  — *  j  transition,  and  a  is  the  upconversion  coefficient, 
which  has  dimensions  of  volume.  I  is  the  intensity  of 
incident  light,  and  Tc  is  the  cavity  lifetime. 

To  simulate  the  dynamics  in  a  cavity  after  illumination 
stops,  the  stationary-state  populations  reached  in  levels 
|0)  to  |3>  during  constant  pumping  were  obtained.  Then 
we  solved  the  undriven  equations  for  times  after  pumping 
ceased  by  utilizing  the  stationary-state  parameters  Q(0) 
and  iVi(O)  as  initial  conditions.  The  resulting  level  \2) 
population  dynamics  are  shown  in  Pig.  3. 

Two  aspects  of  the  dynamics  calculated  with  this 
rate  equation  approach  are  interesting.  First,  for  times 
shorter  than  the  lower-state  lifetime  but  greatly  in  excess 
of  the  cavity  lifetime,  it  is  possible  for  a  population  inver¬ 
sion  to  persist  between  levels  |2>  and  ll>  after  pumping 


has  stopped.  When  the  net  rate  of  pumping  into  state 
|1)  drops  to  zero,  the  rate  of  cooperative  upconversion 
that  feeds  the  upper  laser  level  remains  nearly  constant 
initially,  so  the  inversion  improves  for  a  short  time.  Sec¬ 
ond,  the  optical  feedback  in  a  cavity  introduces  transient 
atom -field  interactions  that  cannot  otherwise  take  place 
(see  Subsection  2.C).  Oscillations,  reflecting  population 
pulsations  in  the  upconversion  medium,  are  predicted  in 
the  rate  equation  limit  and  should  be  evident  in  light 
emitted  from  the  cavity. 

In  the  absence  of  inversion  or  a  cavity,  such  population 
pulsations  are  possible  only  in  spatially  coherent  coopera¬ 
tive  interactions.  To  verify  this^  the  role  of  delocaliza¬ 
tion  in  excited  states  must  be  examined  carefully  with  a 


Time  (ms) 

(a) 


Time  (ms) 
(b) 


Fig.  3.  Numerical  rate  equation  calculation  of  pulsations 
in  the  inversion  of  the  paii^pumped  laser  after  termination 
of  steady-state  excitation.  There  are  no  free  parameters: 
intensity  I  *  153W/cm*  (1.7  times  threshold)  in  a  cav¬ 
ity  (tc  *  1.5  ns)  of  dopant  density  6.7  x  10*^  cm"*  with 
y%  »  0.0714  m8“^,  ys  *  0.179  ms"'^,  ysi  *  0.0780  ms"^, 
732  =“  0.103  ms'S  74  -  143  ms”^,  Boi  =*6.9  x  10 cmVJ, 
Si2  =  2.4  X  10^^®  cmVsr^®  and  a  *  3.2  x  cm*  (a) 

The  initial  part  of  postexcitation  decay  ,  (inset,  magnified  trace 
with  temporal  resolution  on  the  microsecond  time  scale  thowing 
quasi-periodicity),  (b)  Postexcitation  decay  at  coarse  temporal 
resolution  [inset,  fast-Fourier-transform  (FFT)  spectrum  1. 
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Atom  A  Atom  B  Pair  Statas  |A>|B> 

4.  Energy-level  diagram  of  a  model  pair  ayatem.  Left, 
individual  atoma  A  and  B  may  be  pictured  aa  each  abaorbing  one 
pump  photon  (thin  vertical  arrow)  to  reach  the  interaction  level 
indicated  by  the  fiUed  circlea.  Pair  upconversion  then  leada 
to  upconversion  fluorescence  after  accommodation  of  the  energy 
mismatch  A,  aa  indicated  by  the  thick  vertical  arrow.  Right, 
an  alternative  representation  ahowa  the  coupled  pair  system 
aa  a  ladder  of  product  states  |A)|B).  Absorption  d  two  pump 
photons  is  required  for  promotion  of  the  pair  system  to  state  111). 
Direct  excitation  to  state  |02)  by  two-photon  absorption  of  a  sin^e 
atom  ia  off  resonant  because  of  the  energy  defect  A.  Emiaaion 
from  state  |02)  occurs  after  spontaneous  pair  upconversion,  which 
may  be  viewed  as  an  internal  relaxation  from  state  111)  to  |02) 
accompanied  by  emission  of  phonons  for  energy  conservation.. 

different  formalism.  Althonifo  energy  transfer  is  widely 
treated  with  rate  equationa  as  a  one-way,  incoherent 
process  between  a  donor  and  an  acceptor,’  a  density- 
matrix  treatment  can  be  used  to  incorporate  the  possibU- 
ity  of  coherent  back  transfer.  Following  this  approach, 
we  investigated  the  combined  consequences  of  optical  and 
spatial  coherences  in  upconversion  within  the  interacting 
pair  system  shown  schematically  in  Fig.  4.  Analytic  re¬ 
sults  showed  that  sustained  pulsations  are  expected  in 
the  populations  of  excited  states  in  the  intermediate  to 
strongly  coupled  regime  as  the  result  of  coherent  coopera¬ 
tive  upconversion. 

Cooperative  upconversion  results  from  the  exdted-state 
atom- atom  interaction  between  (unaymmetrized)  product 
states  |02)  and  111)  in  the  pair  manifold  shown  in  Fig.  4. 
Interactions  between  other  excited  pair  states  such  as 
|01)  and  |10>  are  prevalent  in  rare-earth  systmns  too  but 
give  rise  only  to  companion  phenomena  socb  as  rapid 
energy  migration  between  sites.  Such  pgacesses  do  not 
contribute  to  upconversion  and  are  neg^acthd  hare.  Also, 
we  considered  the  systen  at  first  to  be  imdrhren.  Inci¬ 
dent  light  prepares  the  qwto  at  timo  f  0  in  the  doubly 
excited  pair  state  111)  M^has  zero  intensity  thereafter. 
Modifications  of  energy-transfer  processes  in  the  presence 
of  radiation  are  considered  in  Subsection  2  J). 

Exdted-state  djmamica  may  bh  calculated  fttan  the 
Ldouville  equation  for  the  density  matrix  p,  augmented 
by  relaxation  terms: 

^ iftyp.  (9) 

The  Hamiltonian  ia  taken  to  he  H  ^  +  H^,  where 

Ho  ia  the  iadated  ion  Hamiltonian  and  Hu  i*  the  in* 
teraction  between  two  coupled  atoms.  Typically,  Hu 
arises  from  multipole- multipole  or  exchange  interactions 
between  electrons  on  adjacent  atoms,^’  but  its  detailed 
form  is  not  of  concern  from  the  point  of  view  of  basic 


dynamics.  Relaxation  terms  represented  globally  by  y 
indude  population  decay  of  each  level  i  at  the  rate  yt 
and  dephasing  of  polarization  on  the  i  j  transition  at  the 
rate  yij  «  (y,-  +  r>)/2.  Off-diagonal  elements  of  H  wee 
written  as  state-specific  coupling  parameters  multiplied 
by  Planck’s  constant.  This  permits  the  strength  of  the 
ion-ion  interaction,  for  instance,  to  be  given  as  a  charac¬ 
teristic  firequency  ihL  *  <02|Hintlll),  thus  avoiding  the 
need  to  specify  multipole -multipole  interactions  or  ac¬ 
companying  phonon  processes  explidtly.  When  the  en¬ 
ergies  of  the  initially  uncoupled  states  111)  and  |02)  are 
degenerate,  one  finds  their  eigenvalues  in  the  presence  of 
interaction  by  diagonalizing  the  subdeterminant  of  matrix 
elements  of  H.  For  example,  with  the  pair  basis  shown 
in  Fig.  4,  the  eigenvalues  of  suhmatrix 

„  [  E  l 


are  readily  found  in  the  subspace  of  states  111)  and  |02) 
to  be 

(11) 

revealing  an  energy  splitting  of  the  doubly  exdted  eigen¬ 
states  equal  to  2h(LL*y^.  This  is  the  well-known 
Davydov  splitting.  The  corresponding  eigenstates  of 
the  interacting  pair  are  4^ 

We  examined  the  temporal  development  of  superposi¬ 
tion  states  rather  than  stationaiy  states.  States  suefo  as 
|00),  111),  and  |02)  are  spatially  localized,  and  their  dy¬ 
namics  directly  reveal  the  spatial  flow  of  energy  within 
the  system.  This  feature  is  extremely  useful  for  imder- 
standing  and  predicting  the  energetics  of  the  system.  Be¬ 
cause  these  states  furnish  a  complete,  alternative  basis 
for  describing  pair  dynamics  in  which  intuitive  ideas  are 
easily  incorporated,  all  the  calculations  were  performed 
with  them. 

Introducing  the  convenient  notation  |00)  »  |0),  111)  » 
|1),  and  |02)  —  |2),  one  can  immediately  write  equations 
of  motion  for  the  pair  superposition  states  in  the  absence 
of  driving  fields  as  follows.  Assuming  the  Davydov  split¬ 
ting  ia  energetically  negligible,  one  finds  that 

^  yiPn  f  (13) 


'  (Lpn  +  E*pu)  -  y%pt$ , 


-  Upn  -  pti)  ^  niPi*  A 
The  anatsrtie,  steody-^te  iidivtuiai  are  .' 


..  ;  .  M  ?  r  .fe  . 

Pn(t)  - 

Pait)  »  exp(-yur)  j[A  (1  -  cos  uat] 


COS  ead 


pia(f)  -  exp(-yia/) 


a>o  Bin 


J[A  (1  -  cos  Wot) 


(18) 
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In  these  expressions  we  have  used  the  following  defini¬ 
tions: 

A  =  (ri-r2)/2,  (19) 

a>o  =  (4LZ-*- (20) 
Pl2^(0  =  Lp2l  +  L* p  12  •  (21) 

Figure  5  shows  the  general  behavior  of  these  solutions 
for  assumed  initial  conditions  pu  ~  1,  p22  =  0,  and  P12'  = 
0.  The  values  of  key  parameters  approximate  those  of 
Tm^"^  in  the  crystal  host  YAIO3,  for  which  relevant  re¬ 
laxation  rates  are  available  in  the  literature.^®  Two  val¬ 
ues  of  the  unknown  coupling  strength  L  were  considered. 
For  L  =  3  kHz  (three  times  the  assumed  dephasing  rate), 
pulsations  in  the  populations  of  pair  superposition  states 
|1)  and  1 2)  appeared  as  sinusoidal  variations  of  pn  and 
P22  versus  time.  For  L  =  1  kHz  (equal  to  the  dephasing 
rate),  no  oscillations  in  population  were  predicted. 

The  population  pulsations  shown  in  Fig.  5(a)  arise 
firom  periodic  reversal  of  the  cooperative  upconversion 
process,  a  manifestation  of  (spatially  coherent)  partial 
delocalization  of  the  excitation.  Theoretically,  the  oscil¬ 
lation  frequency  provides  a  direct  measure  of  the  cou¬ 
pling  strength  L  through  the  relation  for  o)q  given  in 
Eq.  (20).  At  present  this  parameter  may  be  estimated 
only  roughly  fi:x)m  measured  pair  cross-relaxation  rates. 
However,  this  procedure  does  yield  a  value  L  =  1  kHz, 
which  coincides  with  the  experimental  dephasing  rate 
of  level  |1)  (corresponding  to  the  ^^4  level  of  Tm^"^)  in 
Tm:YA103.^®  Because  the  coupling  and  decay  rates  are 
the  same,  this  estimate  suggests  that,  experimentally,  L 
may  be  at  the  critical  value  separating  regimes  of  local¬ 
ized  firom  coherently  delocalized  dynamics.  In  Tm:YA103 
or  other  rare-earth  crystals  effects  of  delocalization  may 
therefore  be  observable  without  optical  cavities.  Weakly 
delocalized  dynamics  may  even  seed  unstable  behavior 
in  the  presence  of  gain,  an  idea  explored  in  more  detail 
in  Subsection  2.C. 

C.  Instabilities  in  Cooperative  Djnnamics 
Because  upconversion  is  inherently  a  nonlinear  process, 
it  can  in  principle  exhibit  unstable  behavior.^  Instabili¬ 
ties  in  excited-state  populations  or  chaos  might  therefore 
be  expected  to  arise  from  interatomic  coupling  of  atoms^^ 
in  at  least  two  ways.  First,  one  can  imagine  cooperative 
upconversion  initiating  population  relaxation  osi^ations 
of  the  type  indicated  in  FHg.  3,  which  might  be  amplified  or 
sustained  in  the  steady  state  by  a  driving  field  in  a  cavity. 
In  this  case  oscillations  would  arise  fiom  cavity- atom  in¬ 
teractions  accompanied  by  stimulated  emission.  Second, 
the  delocalization  disoissed  in  Subsection  2.B  could  initi¬ 
ate  population  pulsations  and  might  be  strong  enough  to 
mediate  an  instability  even  without  a  cavity.  In  this  case 
unstable  behavior  would  arise  purely  firom  atom— atom 
interactions  driven  by  an  optical  field,  without  stimulated 
emission.  In  this  subsection  we  find  the  surprising  re¬ 
sult  that  coherent  delocalization  is  necessary  to  sustain 
steady-state  oscillation  with  or  without  a  cavity. 

As  an  initial  step,  imstable  solutions  of  the  rate  equa¬ 
tions  (4)-(8)  were  sought.  Three-level  lasers  with  lin¬ 
ear  driving  terms  respond  universally  to  departures  from 
steady-state  intensity  by  exponential  decay  back  to  the 


original  steady  state.^  However,  upconversion  is  a  non¬ 
linear  process,  and  a  reexamination  of  this  result  in  a 
multilevel  pair  system  is  warranted. 

A  linear  stability  analysis  of  Eqs.  (4)-(8)  was  per¬ 
formed  to  investigate  cooperative  upconversion  in  a  cavity 
in  the  limit  of  weak  interatomic  coupling.  The  steady- 
state  population  (pii)„  for  each  excited  state  |i)  was 
found  for  fixed  intensity.  Then  linear  departures  from 
steady-state  values  were  substituted  into  the  dynamical 
equations  through  the  replacement  pu  —  (p^j)„  +  Apjj. 
Eigenvalues  of  the  equations  of  motion  of  Apjj  were  then 
examined  to  determine  whether  deviations  fix)m  steady- 
state  values  grew  or  decayed  with  time.  This  analysis 
(see  Appendix  A)  revealed  that  populations  always  relax 
back  to  their  steady-state  values,  with  or  without  a  cav¬ 
ity,  for  arbitrary  values  of  incident  intensity  and  decay 
constants.  That  is,  rate  equation  analysis  leads  to  the 
conclusion  that  pair-pumped  upconversion  lasers  are  uni¬ 
versally  stable. 

In  the  intermediate-to-strong  coupling  regime  a  simi¬ 
lar  stability  analysis  of  upconversion  dynamics  requires  a 
density-matrix  calculation.  This  begins  with  augmenta¬ 
tion  of  the  unperturbed  Hamiltonian  by  the  field  Hamil¬ 
tonian.  In  the  basis  of  uncoupled  pair  states  |0),  11),  and 
1 2)  this  yields 

*  £^|1)<1|  +  E|2)(2|  +  \a\^h(a\a){a\ .  (22) 

Here  \a)  is  a  coherent  state  of  the  radiation  field,  chosen 
to  represent  the  atomic  interaction  with  a  laser  field. 


0  1  2  3  4  5  6 


Time  (ms) 

Fig.  5.  Density-matrix  calculation  of  exdted  pair  state  popula¬ 
tions  versus  time  in  the  presence  of  delocalization.  The  initial 
state  is  characterized  by  pu  «  1,  p22  *  0,  n  «  >2  "  1  itHz, 
and  two  values  of  interatomic  coupling  are  considered,  namely, 
(a)  L  «  3  kHz  and  (b)  L  ^  1  kHz.  Population  pulsations  are 
evident  for  the  larger  value  of  coupling. 
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and  lap  =  (a|a^a|a)  is  the  expectation  value  of  the  num¬ 
ber  operator.^  The  energy  contribution  from  vacuum 
fluctuations  has  been  dropped.  The  ion- ion  interaction 
Hamiltonian  in  this  subspace  is  augmented  by  fleld- 
induced  dipole  transitions  between  the  ground  state  |0) 
and  the  doubly  excited  pair  state  |1)  (not  state  |2»: 

=  iftL|l)<2|  +  h.c.  +  ID (01  +  h.c.  (23) 

In  this  expression  Ato  is  the  |00)  — ►  |01>  transition  fre¬ 
quency,  and  h.c.  denotes  the  Hermitian  conjugate.  The 
parameter  f  is  the  dipole  moment  on  the  |0)  — ►  |1>. transi¬ 
tion,  defined  by  fa  =  {pL'  E}/h,  Only  the  double-atom- 
excitation  terms  of  interest  have  been  retained  from  the 
second-order  Hamiltonian  +  Ma)  *  E101)(01|(/aa  + 

fjLa)  *  E\00)/hAa).  The  subscripts  A  and  B  refer  to  atoms 
A  and  B  of  a  coupled  pair.  Two-photon  excitations  of 
individual  atoms  have  been  specifically  excluded  from 
Eq.  (23). 

The  populations  of  the  pair  system  in  the  interaction 
representation  are 


dpoo 

dt 

dpn 

dt 

dp22 

dt 


==  P3  +  ripii  +  r2P22 , 


P3  +  P4  -  riPii , 


=  -P4  -  72P22  * 


(24) 


(25) 


(26) 


The  equations  for  optical  and  spatial  coherences  in  the 
system  are  given  by 


dp3  ^ 
d^ 


i\pn  ~~  Poo)  “  Toip3 » 


dp4 

dt 


=  2y3(p22  pii)  ^  ri2p4, 


(27) 


(28) 


respectively.  Here  yoi  and  yx2  are  the  optical-  and  the 
spatial-coherence  dephasing  rates,  and  we  have  used  the 
following  definitions: 


P3  =  (/^«)Vo  +  (r«Vpoi, 

P4  =  -L(pi2  P2i), 


(29) 

(30) 

(31) 

(32) 


We  emphasize  the  role  of  incident  li^t  by  replacing  the 
expectation  value  of  the  number  operator  (a|n|a)  =  ta|* 
by  a  quantity  I  proportional  to  (but  not  identical  with) 
intensity,  as  defined  in  £q.  (31). 

In  these  equations  the  coherence  term  po2,  which  is  sec¬ 
ond  order  in  the  field,  has  been  ne^ected.  This  at  first 
appears  ui^ustified  because  the  coherence  term  poi,  which 
is  similarly  second  order  in  the  pumping  field,  has  been 
retained.  However,  it  can  easily  happen  that  sequential 
or  simultaneous  two-photon  absorption  processes  fail  to 
populate  doubly  excited  states  while  cooperative  upcon- 
version  takes  place.  If  phonons  accommodate  the  energy 
defect  A,  upconversion  can  occur  rapidly  and  may  gener¬ 
ate  weak  spatial  coherence  or  the  delocalization  described 


by  poi*  At  the  same  time,  if  energy  levels  are  not  in 
exact  coincidence  with  the  incident  frequency,  two-photon 
optical  coherence  does  not  develop  ( po2  =  0). 

In  a  closed  three-level  pair  system,  population  is  con¬ 
served.  Hence  we  can  eliminate  groimd-state  variables, 
using  the  relation  poo  +  Pii  +  P22  “  1-  Making  the  sub¬ 
stitutions  pii  — ►  Pi  and  P22  P2  for  convenience,  we  find 
that  the  equations  of  motion  in  matrix  form  become 


Pi 

P2 

PS 

>4  J 


-n  0  1  -1 

0  -72  0  ”1 

2/"  "701  0 

l-2p  2P  0  -7uJ 


Pi 

P2 

PS 

Pa 


*0' 

0 

/2 

0 


(33) 


The  overdots  indicate  differentiation  with  respect  to  time. 
The  secular  equation,  of  the  general  form  +  osA^  + 
azX^  +  uiA  +  oo  =  0,  has  particularly  simple  coefficients 
if  level  11)  is  assumed  to  be  metestable  (yi  »  0).  Then 
the  secular  determinant  factors  into  a  trivial  solution  with 
zero  eigenvalue  and  a  cubic  equation  that  can  be  solved 
analytically.  Here  we  are  interested  only  in  identifying 
a  condition  that  renders  the  real  part  of  any  eigenvalue 
positive  because  then  the  system  will  oscillate  contin¬ 
uously.  The  Routh-Hurwitz  criterion^  for  such  a% 
instability  is  that  the  system  determinant  D  —  0102^3  <¥ 
ai2  "  ooosa  *  i572*(16^  +  “  3fl/2  must  be  zero  o§ 

negative.  This  condition  occurs  when  the  intensity  | 
exceeds  the  critical  value 

3 


(34) 


From  this  result  it  may  be  concluded  that  spatial  coher¬ 
ence  can  generate  population  oscillations  in  the  steady 
state  without  optical  feedback.  This  feature  of  coherent 
cooperative  upconversion  has  not,  to  our  knowledge,  been 
previously  discussed  and  may  be  a  useful  new  signature  of 
delocalization  in  upconversion.  A  numerical  calculation 
of  the  growth  of  this  instalnlity  is  shown  in  Fig.  6. 

D.  Enhanced  Upconversion  and  Energy  Transfer 
We  now  show  that  the  energy-trazi^r  rate  between 
atoms  is  frequency  dependmt  when  modulated,  incoher¬ 
ent  light  is  used  to  excite  pair  lipoonversiorL  A  resonant 
condition  for  enhanced  upconversion  emerges  from  this 
analysis. 

The  dynamical  equations  that  we  consider  are  identical 
to  Eqs.  (13)-(15),  except  that  an  incoherent  driving  term 
is  added  to  the  ri^it-hand  side  of  Eq.  (14).  ’  The  modified 

(35) 

•a#-,:..'.  - 

The  optMxt^piraiiiuvraie  /i^  de^  inten- 

sily,  the  Bnetafa  ggcbe<WHWf:iB#^iep<yulatio^  dif. 
ference  betweefl  lente  |^«nt  ^^'ltt-^le;iow.<lepletion 
limit  of 

the  dauil^NBati^ '<^mMtai^0WlMigic1^  equations  of 
motion  in  toms  of dB<fiio4'h»  B^(2IB^  'ne  find  that 


(36) 


Pu  ” 

’“71 

0 

1  " 

"Pii  ’ 

"  0  “ 

P22 

= 

0 

-ra 

-I 

Paa. 

Ap(t) 

.P12'- 

.-2/8 

“7U- 

.Pia'_ 

_  0  _ 

P.  Xie  and  S.  Rand 


Vol,  11,  No.  5/May  1994/J.  Opt.  Soc.  Am.  B  907 


^  0.48 


*<s 

eQ 

3 

a 

o 

•o 

V 

N 


0.46  r 


^  0.44 

eg 

E 


o 

Z 


0.42  h 


xlO- 


4  6 

Time  (ms) 


Fig.  6.  Growth  of  quasi-periodic  output  from  a  constant 
steady-state  condition  of  a  pair-pumped  laser  (density-matrix 
calculation),  (a)  —  2.6  x  l(r  kHz^,  ^  =  900  kHz^,  and 

72  —  100  cm~^,  71  —  0.  (b)  Fourier  spectnim  of  (a). 

If  the  incident  field  is  modulated  at  a  frequency  o)  ac¬ 
cording  to  Ap(i)  =  Ap(0)exp(i<y0,  then  harmonic  solutions 
exist  of  the  form 


~  Pllit)  ~ 

rpu(o)i 

P22(f) 

P22(0) 

_Pl2'(f)  _ 

^12'(0)_ 

exp(ia;^) , 


(37) 


and  substitution  into  Eq.  (36)  yields  the  solutions  given  by 


P22(0) 


2P 


(•yi2  +  iai)\_{yi2  +  i<u)*  +  oi»o^] 


Ap(0),  (39) 


_ -2/3(72  +  i(») _  . 

(yi2  +  »<«)[(ri2  +  ’’ 


(40) 


In  these  expressions  a  resonant  frequency  too  appears, 
defined  by 


too  — 


4/3 


1/2 


(41) 


The  physical  meaning  of  Eqs.  (38)-(40)  can  best  be  ap¬ 
preciated  by  comparison  of  the  fluorescence  intensity  at  a 
particular  fi-equency  with  that  at  zero  firequency  (a>  =«  0), 
with  a  response  function  defined  by 


R(w)^ 


P22(0) 

P22(0)L-0 


712(712^  +  ^0^) 


{(iWO^  +  712^)[{<Oq^  -  <y2  +  +  ^f^2y^^2YJ2 


(42) 


R{<m})  is  proportional  to  emission  intensity,  which  is 
dependent  on  energy  transfer,  and  is  plotted  in  Fig.  7  for 
two  values  of  the  interatomic  coupling  strength  L,  again 
with  the  use  of  decay  parameters  typical  of  Tm^YAlOa. 
For  weak  coupling  (L  =  0.5  kHz),  the  energy-transfer  re¬ 
sponse  function  falls  off  rapidly  as  the  frequency  of  modu¬ 
lation  is  increased  from  zero.  However,  for  stronger 
coupling  (L  —  5  kHz)  it  exhibits  a  resonant  peak  near  the 
coupling  frequency.  The  precise  frequency  at  which  the 
peak  occurs  is  given  by 


^2  _ 


-  4yi2^  +  -  32010^712^)^ 


(43) 


When  the  dephasing  rate  of  spatial  coherence  is  much 
less  than  the  coupling  firequency,  the  energy-transfer  reso¬ 
nance  is  well  de^ed.  When  the  dephasing  rate  exceeds 
the  coupling  fi^uency,  the  peak  broadens  and  becomes 
less  pronounced,  eventually  disappearing. 

3.  EXPERIMENTAL  PROCEDURES 

To  investigate  some  of  these  predictions,  experiments 
were  performed  in  single  crystals  of  5%  Er:CaF2  and  5% 
Er:LiYF4.  Two  sample  disks  of  3-mm  thickness  were 
prepared  as  monoUtMc,  hemispherical  optical  cavities, 
and  a  third  sample  of  the  same  thickness  was  polished 
plane  parallel.  The  monolithic  crystals  were  cut  and 
were  polished  identically  for  laser  experiments,  with  one 
flat  surface  and  one  convex  surface  of  radius  2.5  cm. 
Both  surfaces  were  antireflection  coated  in  the  range 
1.4- 1.6  /zm.  Most  observations  were  made  in  the  first 


Fig.  7.  Response  function  for  enhanced  energy  transfer  result¬ 
ing  from  modulation  of  the  optical  driving  field  for  two  values 
of  the  interatomic  coupling  strength  L,  (a)  L  »  5  kHz,  (b) 
L  *  0.5  kHz.  Other  parameters  used  are  71  —  0.43  kHz  (twice 
the  measured^®  decay  rate),  72  *  1.6  kHz,  and  712  =  1-02  kHz. 


908  J.  Opt.  Soc.  Am.  B/Vol.  11,  No.  5/May  1994 


P.  Xie  and  S.  Rand 


crystal,  in  which  the  curved  surface  coating  gave  total 
reflection  between  2.7  and  2.9  fim  and  the  flat  served 
as  a  2%  output  coupler.  The  other  monolithic  sample 
was  coated  for  output  at  854  nm.  The  third  sample  was 
used  in  an  astigmatically  compensated,  three-mirror  cav¬ 
ity  designed  for  operation  at  702  nm.  Collectively  these 
samples  furnished  light  sources  in  which  emission  was 
sustained  by  pair,  trio,  and  quartet  upconversion  pro¬ 
cesses,  respectively. 

Each  sample  was  pumped  longitudinally  with  a  cw 
NaCl  color-center  laser,  absorbing  approximately  75% 
of  incident  radiation  at  1.51  ^m.  Laser  action  at 
2.8  fim  was  achieved  at  room  temperature,  as  described 
previously.  The  lasers  with  output  wavelengths  at  855 
(Ref.  25)  and  702  nm  (Ref.  26)  required  cooling  to  liquid- 
nitrogen  and  liquid-helium  temperatures,  respectively. 
Cw  laser  emission  was  observed  in  each  sample  with  out¬ 
put  powers  in  the  milliwatt  range,  sufficient  for  a  variety 
of  experiments  related  to  cooperative  dynamics.  Laser 
output  of  the  first  crystal  was  monitored  through  a  notch 
filter  at  2.8  pum  (FWHM  10  nm)  with  a  fast  InAs  photo¬ 
diode  (rise  time  —5  ns).  Dynamics  at  854  and  702  nm 
were  monitored  with  fast  photomultipliers  and  were 
recorded  on  a  l-GHz  digitid  oscilloscope.  Thermopile 
detection  was  used  for  power  measurements. 

Two  main  types  of  time-resolved  measurement  perti¬ 
nent  to  a  comparison  with  the  theoretical  concept  out¬ 
lined  in  Section  2  were  made.  The  first  consisted  of 
observations  of  2.8- Aim  laser  output  transients  at  the  lead¬ 
ing  and  trailing  edges  of  square  pump  pulses  generated 
by  acousto-optic  techniques.  We  made  single-shot  digital 
recordings  with  100-ns  temporal  resolution  (analog  band¬ 
width),  using  a  1-MHz  sampling  rate  over  8,000  points. 
Second,  output  intensity  was  digitally  recorded  versus 
time  in  the  instability  regime,  with  constant  pumping. 


these  laser  measurements  whether  cooperative  enhance¬ 
ment  affected  performance. 

Nevertheless,  other  recent  observations  support  the 
concept  of  cooperatively  enhanced  quantum  efficiency  out¬ 
lined  in  Subsection  2.A.  Output  power  measurements  of 
a  2.8- Aim  laser  pumped  by  conventional  means  at  980  nm 
in  concentrated  Er:LiYF4  recently  revealed  performance 
marginally  in  excess  of  the  theoretical  efficiency  limit 
set  by  rfE  =  Am/Aout  —  35%.^^  This  result  was  ascribed 
to  pair  upconversion  dynamics,  furnishing  the  first  ex¬ 
ample  of  a  laser  with  output  enhanced  by  cooperative 
dynamics.  Enhancement  should  extend  to  cooperative 
upconversion  lasers  as  well,  however.  Although  our  own 
measurements  do  not  reveal  enhancement  of  upconver¬ 
sion  laser  output,  we  do  furnish  direct  evidence  of  recy¬ 
cling  dynamics  in  Subsection  4.B. 


Pump  Power  (mW) 
(a) 


4.  RESULTS  AND  DISCUSSION 


A.  Enhanced  Quantum  Efficiency 
The  most  efficient  upconversion  laser  in  our  experi¬ 
ments  was  the  imoptimized  trio  laser,  which  operated 
at  855  fim.  Measurements  yielded  a  slope  efficiency 
of  28%  and  an  overall  energy  efficiency  of  26%  for 
this  laser.  The  latter  value  is  to  be  compared  with 
the  maximum  theoretical  limits  with  and  without  re¬ 
cycling,  which  give  rfg*  —  [l/(m  -  l)]Ain/Aout  “  88%  and 
ris  ^  (l/m)Xin/ X^i  «  59%,  respectively.  Although  the  cw 
trio  laser  worked  remarkably  well  on  a  self-terminating 
transition,  the  actual  laser  efficiency  gave  no  indication 
of  quantum-efficiency  enhancement. 

Figure  8  shows  representative  data  for  upconversion 
laser  emission  on  the  transition  at  701.5  nm, 

excited  by  1.5- Atm  radiation  and  originating  from  a  state 
at  nearly  four  times  the  incident  photon  energy.  The 
mere  operation  of  this  fourfold  upconversion  laser  is  sur¬ 
prising  and  immediately  raises  questimis  as  to  what  fun¬ 
damental  limits  may  exist  for  cw  laser  action  involving 
high  degrees  of  upconversion  and  whether  enhancement 
mechanisms  are  operative  in  such  devices.  However,  the 
fourfold  upconversion  laser  was  also  less  efficient  than  ei¬ 
ther  the  maximum  value  tje*  71%  or  tie  ^  53%.  Con¬ 
sequently  it  was  not  possible  to  conclude  directly  from 
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Fig.  8.  Output  versus  input  power  for  the  fourfold  cw  upconver¬ 
sion  laser  operating  at  701.5  nm  (A«x  1>5  Atm).  Inset,  laser 
emission  spectrum,  (b)  Energy  diagram  identifying  the 
levels  involved  in  the  pumping  and  the  emission  processes  of 
the  fourfold  upconversion  laser. 
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(b) 

Pig.  9.  Oscillations  in  output  power  observed  at  (a)  the  leading 
and  (b)  the  trailing  edges  of  a  2.8-/Am  Er^LiYF4  laser  pulse 
excited  by  a  square  pulse  at  =  1.5  fim  of  10-ms  duration. 
Excitation  was  gated  acousto-optically  from  the  output  of  a  cw 
NaCl  laser.  The  rise  and  fall  times  of  the  pump  pulses  were 
less  than  50  ns.  Insets,  Fourier  transforms  showing  dominant 
frequency  components. 


B.  Population  Pulsations 

When  cw  operation  of  the  2A^fim  pair-pumped  laser  was 
terminated,  a  transient  increase  in  output  power  was  ob¬ 
served.  This  result  is  in  accord  with  the  simple  rate 
equation  treatment  presented  in  Subsection  2.B.  In  ad¬ 
dition,  for  excitation  at  A«z  »  1.5  /tim,  unexpected  rapid 
oscillations  occurred  under  the  envelope  of  thin  tran¬ 
sient  peak,  as  shown  in  Fig.  9.  Relaxation  oscillations 
were  observed  at  the  leading  edge  of  pump  pulses  also 
[Fig.  9(a)]  but  are  not  discussed  here  because  of  the  ana¬ 
lytic  complication  of  including  a  driving  field  in  the  dy¬ 
namics.  Conclusions  may  be  drawn  more  straightfor¬ 
wardly  firom  an  examination  of  trailing-edge  transients. 

An  enlargement  of  postexdtation  pulsations  is  shown 
in  Fig.  9(b),  together  with  their  Fourier  spectrum,  reveal¬ 
ing  a  chirped  frequency  spectrum  centered  near  210  kHz. 
This  frequency  is  very  near  the  measured  inverse  lifetime 
(I/ttt)  -  210  kHz  of  the  V9/3  level  of  Er3'':LiYF4.  From 
Pig.  9(b)  it  is  clear  that  pulsations  of  the  stimulated  emis¬ 


sion  continued  long  after  external  excitation  of  the  system 
ceased.  This  clearly  indicates  that  they  originate  from 
the  interaction  of  the  residual  atomic  inversion,  the  op¬ 
tical  cavity,  and  the  cavity  photon  density.  Moreover, 
these  pulsations  are  in  excellent  quantitative  agreement 
with  the  rate  equation  calculation  given  in  Fig.  3,  with 
no  adjustable  parameters. 

The  origin  of  these  transient  pulsations  can  be  pictured 
in  the  following  way.  Once  external  excitation  ceases, 
inversion  can  recur  only  through  pair  upconversion  to 
the  ^/9/2  level  (analogous  to  level  |3>  in  Fig.  2)  of  Er^^. 
This  process  is  followed  by  nonradiative  relaxation  to 
the  upper  laser  level  (|2>  in  Fig.  2).  While  the  up¬ 
conversion  rate  is  very  high,  the  limiting  rate  at  which 
inversion  can  be  reestablished  following  depletion  of  the 
gain  is  set  by  the  inverse  lifetime.  As  long  as  the 
level  (|1)  in  Fig.  2)  responsible  for  pair  upconversion 
remains  well  populated,  cooperative  upconversion  efiB- 
ciently  replenishes  the  upper-state  population  in  this  way. 
Inversion  density  builds  until  the  threshold  is  reached, 
whereupon  stimulated  emission  depletes  the  accumulated 
gain.  Quite  distinct  firom  conventional  relaxation  oscil¬ 
lations,  these  popvdation  pulsations  are  a  direct  manifes¬ 
tation  of  population  recycling  by  cooperative  dynamics, 
amplified  by  the  feedback  within  an  optical  cavity. 

C  Instabilities  in  Driven  Systems 
Whereas  true  cw  operation  of  the  pair-pumped  2.8-/Urm 
laser  was  achieved  at  low  pumping  intensity,  instabilities 
became  evident  at  intensities  well  above  the  threshold  for 
cw  lasing.  These  results  are  shown  in  Figs.  10-12. 

For  an  excitation  wavelength  of  **  1.5  /^m,  inversion 
at  2.8  fim  is  sustained  purely  by  .cooperative  dynamics.^* 
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Fig.  10.  Traces  of  pair-pumped  Er^CaFa  laser  output  at  2.8  /im 
at  various  pumping  intensities,  (a)  The  pump  pulse,  (b) 
Threshold  operation,  (c)  True  cw  operation  with  leading- 
and  trailing-edge  transients  (1.1  times  threshold),  (d)  Cw 
operation  at  higher  power  (1.5  times  threshold),  (e)  The  growth 
of  sustained  oscillation  fiom  noise  (two  times  threshold),  (f) 
Unstable,  quasi-periodic  output  (three  times  threshold). 
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Frequency  (kHz) 

Fig.  11.  Fourier  transforms  of  the  Fig.  10  time  series,  showing 
the  presence  of  several  modulation  frequencies  at  high  pump¬ 
ing  intensity  but  no  subharmonics.  Leading-  and  trailing-edge 
transients  were  excluded  from  the  analysis. 


Time(iiis) 


Fig.  12,  Exj^rimental  trace  of  the  growth  of  the  instability 
shown  in  Fig.  10(e)  with  high  temporal  resolution.  Inset, 
Fourier-transform  spectrum. 

At  other  excitation  wavelengths  the  cooperative  process 
occurs  at  a  much  reduced  rate  because  the  steady-state 
^/i3/2  population  from  which  pair  upconversion  originates 
is  greatly  reduced.  The  main  result  of  this  section  is  that 
instabilities  are  observed  only  when  the  excitation  wave¬ 
length  is  chosen  to  maximize  the  cooperative  dynamic  con¬ 
tribution  to  laser  inversion  (A«,  =  1.5  /Am). 

The  experimental  observation  of  spontaneous  growth 
of  an  instability  from  noise  above  a  distinct  threshold 
intensity  is  shown  in  Fig.  10.  The  associated  Fourier 
spectrum  is  presented  in  Fig.  11^  and  a  second  record¬ 
ing  is  shown  in  Fig.  12,  with  higher  temporal  resolu¬ 
tion.  These  experimental  results  are  in  excellent  agree¬ 
ment  with  the  density-matrix  prediction  given  in  Fig,  6, 
At  this  pumping  intensity  laser  output  experiences  a 
hard  loss  of  stability  characteristic  of  a  subcritical  Hopf 
bifurcation,^  suddenly  developing  a  deep,  periodic  modu¬ 
lation  above  the  instability  threshold.  Oscillations  grow 


progressively  more  complex  as  pumping  is  increased,  pos¬ 
sibly  because  of  a  secondary  Hopf  bifurcation  on  a  torus. 
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Fig.  13.  Phase  space  plots  of  unstable  output  from  the  2.8-Aim 
pair-pximped  laser  (a)  6x/dt  versus  <  for  the  xit)  time  sequence 
in  the  3-6-ms  interval  shown  in  Fig.  10(e)  at  two  times  laser 
threshold  (Er:CaF2),  (b)  dx/dt  versus  t  at  three  times  laser 
threshold  (Er-CaF 2),  (c)  dx/dt  versus  t  at  five  times  laser  thresh¬ 
old  (Er:LiYF4). 
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This  progression  is  reflected  in  the  attractor  plots  shown 
in  Fig.  13. 

Although  the  laser  is  renowned  for  being  nois¬ 

ier  than  most  other  solid-state  lasers  under  free-running 
conditions,  to  our  knowledge  no  sustained  oscillations  like 
those  reported  here  have  been  observed  for  excitation  at 
pump  wavelengths  other  than  1.5  fim.  In  our  experi¬ 
ments,  in  addition  to  using  Aex  =  1.5  to  pump  2.8-Atm 
Er:CaF2  and  Er:LiYF4  lasers,  we  tried  alternative  wave¬ 
lengths  of  Aex  —  651-657  nm  and  Aex  =  800  nm.  As  in¬ 
dicated  in  Fig.  14,  the  output  observed  from  the  2.8-Aim 
Er:CaF2  laser  pumped  three  times  over  threshold  at  Aex  = 
652  nm  showed  no  spiking  behavior  during  continuous  op¬ 
eration.  After  termination  of  pumping  at  Ae*  “  652  nm 
this  laser  showed  transient  enhancement  of  laser  out¬ 
put  (an  afterpulse),  but  population  pulsations  beneath  the 
afterpulse  envelope  were  notably  absent.  Similarly,  for 
=  800  nm,  no  sustained  oscillations  were  observed  at 
equivalent  levels  above  threshold,  and  no  postexcitation 
pulsations  occurred. 

At  pump  wavelengths  other  than  1.5  /jLxn,  coop¬ 
erative  dynamics  evidently  served  to  overcome  the 
self-terminating  nature  of  the  2.8-/jtm  emission  but 
contributed  only  weakly  to  upper-laser-level  population. 
Oscillations  in  2.8-AAm  laser  emission  did  not  occur  under 
these  circumstances.  Unstable  operating  conditions 
were  reached  only  when  excitation  was  provided  at 


Er  laser  operating  at  2.8  fim  are  due  to  the  nonlinear 
dynamics  of  cooperative  upconversion  and  can  be  repro¬ 
duced  by  a  density-matrix  theory  of  cooperative  upconver¬ 
sion  that  incorporates  spatial  coherence.  Rate  equations 
successfully  describe  transient  upconversion  behavior 
but  fail  to  reproduce  the  observed  steady-state  instabili¬ 
ties,  thereby  indicating  the  importance  of  spatial  coher¬ 
ence  in  pair  dynamics  of  rare-earth  ions  in  upconversion 
lasers.  Coherent  cooperative  upconversion  can  in  prin¬ 
ciple  mediate  popiilation  pulsations  in  rare-earth  systems 
without  cavities  and  admits  a  measure  of  control  over 
energy  transfer.  These  interesting  prospects  warrant 
further  investigation. 

APPENDIX  A 

To  analyze  the  stability  of  the  pair-pumped  upconver¬ 
sion  laser,  the  system  response  to  small  perturbations  in 
Eqs.  (4)-(8)  must  be  examined.  Solutions  of  the  excited- 
state  populations  that  grow  with  time  are  sought  after  one 
writes  the  time-dependent  population  Niit)  of  level  U>  as 
a  sum  of  its  steady-state  value  Nu  and  its  fluctuation  rii, 
according  to  JV'(^)  =  iV^  /if.  Similarly,  Q(^)  ==  Q,  +  q, 
and  it  is  assumed  that  rii/Nu  «  1  and  q/Q,  «  1. 

Using  no  +  'll  +  +  rta  =  0,  one  obtains  the  following 

reduced  set  of  dynamical  equations  for  the  fluctuations 
themselves: 
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1.5  fjLm,  for  only  then  did  nonlinear  cooperative  dynamics 
dominate  the  inversion  mechanism.  We  take  these  re¬ 
sults,  in  combination  with  the  theoretical  considerations 
of  Subsections  2.B  and  2.C,  to  be  strong  evidence  that 
the  observed  instabilities  of  the  pair-pumped  laser  are 
completely  attributable  to  cooperative  nonlinear  dynam¬ 
ics.  Additionally,  one  may  infer  that  noisy  operation  of 
the  2.8-Aim  Er  laser  pumped  at  wavelengths  other  than 
1.5  fim  originates  from  a  subthreshold  instability  of  the 
cooperative  upconversion  contribution  to  gain  in  the  laser 
me^um.  Improved  operation  of  this  laser  might  there¬ 
fore  be  attainable  with  feedback  control  of  the  ty^  re¬ 
cently  applied  to  stabilize  chaotic  laser  output.^ 

5.  SUMMARY 

Cooperative  nonlinear  dynamics  in  highly  doped  media 
have  important  implications  for  upconversion  and  con¬ 
ventional  solid-state  lasers  alike.  In  this  paper  we  have 
introduced  the  concepts  of  population  recycling  in  coop¬ 
erative  upconversion  and  enhancement  of  quantum  effi¬ 
ciency.  These  concepts  are  relevant  to  an  understanding 
of  the  theoretical  efficiency  limits  of  cw  upconversion 
lasers  and  to  potential  improvements  of  conventional 
lasers  in  dense  media. 

Cooperative  dynamics  give  rise  to  population  pulsa¬ 
tions  and  instabilities  that  can  be  observed  in  optical  cavi¬ 
ties.  We  have  established  that  pulsations  dinring  after¬ 
pulses  and  sustained  oscillations  of  the  cw  pair-pumped 


Solutions  of  Eq.  (Al)  exist  of  the  form 
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Fig.  14.  Output  trace  of  the  2.8- Atm  Er:CaF2  laser  pumped  by 
a  long  square  pulse  in  the  red  spectral  region  (Aex  ^  852  nm). 
At  this  pump  wavelength  output  is  continuous  even  at  three 
times  above  threshold,  and  the  trailing-edge  spike  has  no  rapid 
underlying  oscillations. 


912  J.  Opt.  Soc.  Am.  B/Vol.  11,  No.  5/May  1994 


P.  Xie  and  S.  Rand 


where  A  is  £in  eigenvalue  of  the  matrix  A  of  coef¬ 
ficients  in  Eq.  (Al).  Eigenvalues  are  determined  by 
setting  Det[A  -\U]  =  0,  where  U  is  the  unit  tensor, 
and  this  secular  equation  can  be  written  in  the  form 
04  +  oaA^  +  02  A^  +  OiA  +  oq  =  0.  The  coefficients  Oj 

are  determined  by  explicit  evaluation  of  the  secular  de¬ 
terminant.  The  system  becomes  unstable  if  any  root  is 
positive,  for  then  the  small  initial  perturbation  will  grow 
exponentially  in  time  according  to  Eq,  (A2). 

The  Routh—Hurwitz  criterion^'^  for  stability  requires  an 
examination  of  the  following  determinants: 
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Oo  =  abed  +  yslSabd  +  bed  +  bd(yiQ  +  y2o)]  +  bedyso , 
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02  ”  oc  +  3o6  +  26c  +  2bd  +  20(72  +  73) 

+  6(710  +  720  +  273)  +  c(272  -  2730  +  731  +  2732) 
+  [710(72  +  73)  +  7273], 

03  =  2(a  +  6  +  c)  +  7i  +  72  -f-  73 , 

04  =  1 .  (A4) 


Here  we  have  used  the  definitions  o  =  Boil,  6  —  BuQsr 
e  —  2aNis,  and  d  -  Bi2(iV2»  ~  Nu)  ^  and  we  note 
that  all  these  quantities  and  the  coefficients  given  above 
are  positive  above  laser  threshold.  The  system  will  be¬ 
come  unstable  if  any  determinant  Di(i  =  1,2, 3, 4)  is  less 
than  zero. 

Because  04  =  1  we  note  that  Dz,  so  the  number  of 
determinants  that  must  be  considered  is  reduced  to  three: 


■C^i  ”  oi ,  (A5) 

D2  =0x02  *“  0300 ,  (A6) 

Dz  =  032)2  -  0401^ .  (A7) 

From  Eqs.  (A4)  it  is  evident  that  Oj  >  0  for  all  i.  Hence, 

for  an  instability  to  develop,  it  is  necesssury  that  either 
I>2  <  0  or  Z>3  <  0.  However,  one  can  easily  verify,  by 
expanding  Eqs.  (A6)  and  (A7)  directly  and  by  canceling 
all  negative  terms,  that  D2  and  Dz  are  always  positive 
for  positive  inversion  iN2s-Nu).  Consequently  the  rate 
equations  describing  pair-pumped  laser  operation  are  uni¬ 
versally  stable. 
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We  report  a  new  Er^‘‘‘:LiYF4  cryogenic  upconversion  laser  pumped  by  a  fourfold  upconversion 
process.  Excitation  at  1.5  ^m  results  in  laser  emission  at  701.5  nm  on  a  transition  with  an  upper 
state  at  nearly  four  times  the  pump  photon  energy.  Mechanisms  and  the  concept  of  cooperative 
upconversion  enhancement  are  examined. 


Recent  demonstrations*”^  of  room-temperature  upcon¬ 
version  lasers  have  stirred  interest  in  their  potential  as 
practical  sources  of  short  wavelength  radiation  for  display 
and  data  storage  applications,  as  well  as  for  communica¬ 
tions  and  ultrashort  pulse  generation  at  visible  and  ultra¬ 
violet  wavelengths.^  To  date  however,  little  consideration 
has  been  given  to  the  relative  merits  of  various  fundamen¬ 
tal  mechanisms  for  achieving  upconversion,  basic  limits  to 
achievable  degrees  of  upconversion  (ratio  of  upper  laser 
level  energy  to  incident  photon  energy),  or  to  questions  of 
stability  of  the  highly  nonlinear  pumping  processes  in¬ 
volved  in  these  devices.  In  this  article  we  discuss  two  of 
these  topics  in  the  context  of  new  results  demonstrating  a 
continuous-wave,  fourfold  upconversion  laser.  The  impor¬ 
tant  third  issue  of  stability  of  solid  state  lasers  in  which 
nonlinear  cooperative  dynamics  play  an  important  role  has 
been  considered  by  Xie^  and  will  be  published  separately.® 

Upx:onversion  fluorescence  observed  in  Er;LiYF4  due 
to  irradiation  with  a  continuous- wave  (cw)  NaCl  laser  at 
1.5  /xm  at  liquid-helium  temperature  is  shown  in  Fig.  1  (a). 
The  upconversion  mechanisms  responsible  for  fluorescent 
emissions  near  850  and  550  nm  have  been  studied  in  past 
work  and  arise  from  cooperative  (multi-atom)  energy 
transfer  processes^  when  cw  excitation  is  restricted  to  1.5 
/xm. 

Emissions  near  410,  650,  702  nm  from  higher  lying 
states  have  not  yet  been  studied  as  thoroughly,  but  are 
sufficiently  intense  to  draw  attention  as  potential  laser  can¬ 
didates.  The  two  near-ultraviolet  lines  at  407  and  413  nm 
arise  from  transitions  with  upper  levels  at  roughly  four 
times  the  incident  photon  energy,  as  indicated  in  Fig.  1(b). 
The  fluorescent  emission  at  702  nm  arises  from  Stark  com¬ 
ponents  of  the  ||/2transition. 

Laser  experiments  were  performed  in  a  three-mirror, 
astigmatically  compensated^  cavity  consisting  of  two  5  cm 
radius  total  reflectors  and  an  output  mirror  with  97% 
reflectively  at  702  nm.  A  3-mm-thick  crystal  of  5% 
Er:LiYF4  inserted  at  Brewster’s  angle  within  the  focusing 
arm  of  the  laser  served  as  the  gain  medium.  Its  optic  axis 
was  oriented  parallel  to  the  crystal  surface  in  the  plane  of 
incidence  of  horizontally  polarized  pump  radiation  and  the 
crystal  was  suspended  on  a  cold  finger  in  vacuum.  The 
laser  emission  spectrum  at  701.5  nm,  assigned  to  the 
(l)-^Vi  1/2(3)  transition,*®  and  the  variation  of  output 
power  with  input  are  shown  in  Fig.  2.  The  overall  efficiency 
was  0.06%  and  observed  slope  efficiency  was  0.09%. 

The  excitation  spectrum  of  laser  emission  revealed  a 


one-to-one  correspondence  with  erbium  absorption  wave¬ 
lengths  in  the  1.5  /xm  region.  This  result  is  shown  in  Fig.  3 
and  is  similar  to  earlier  findings  for  pair^  and  trio**  lasers, 
but  contrasts  the  restrictive  wavelength  dependence  ex¬ 
pected  for  multiphoton  upconversion  processes.  It  there¬ 
fore  has  important  implications  for  the  inversion  mecha¬ 
nism  of  the  fourfold  la^.*^ 

For  muhiphoton  absorption  to  be  effective,  several 
ground  and  excited  state  aA)sorption  (ESA)  frequencies 
must  overlap  sharply.  Hence,  excited  state  resonances 
should  appear  in  the  excitation  spectrum  when  this  mech¬ 
anism  is  operative.*^  For  cooperative  upconversion,  only  a 
single  pump  transition  is  relevant,  terminating  anywhere 
within  the  metastable  manifold  in  which  cooperative  inter- 


FIG.  1.  (a)  Upconvmioii  fluorescence  spectrum  observed  in  Er^‘^:LiYF4 
with  excitation  at  Ae,=  1.5  pm  ( rs=  10  K),  uncorrected  for  instrumental 
response,  (b)  Schematic  diagram  of  Er^'*'  energy  levels,  showing  emission 
at  410  and  701.5  nm  from  at  nearly  four  times  the  incident  photon 
energy. 
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FIG.  2.  Output  vs  input  power  of  the  cw  fourfold  upconverston  laser 
operating  at  701.5  nm  (r=sl0  K).  Inset:  Laser  emission  spectrum. 
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FIG.  4.  Mechanism  for  enhanced  quantum  efficiency  in  cooperative  up- 
conversion  systems.  Following  an  initial  excitation  step  (by  direct  pump¬ 
ing  or  upconversion),  a  cascade  from  higher  to  lower  excited  states  occurs 
with  branching  ratios  of  17^  at  each  step.  Wiggly  arrows  indicate  nonradt- 
ative  decay.  A  given  atom  can  emit  more  than  one  photon  per  excita¬ 
tion  if  the  upconversion  rate  exceeds  the  spontaneous  decay  rate  of  the 
lowest  excited  state  (and  detailed  balance  is  satisfied). 


actions  occur.  In  this  case  only  ground  state  absorption 
resonances  should  appc^  in  the  excitation  spectrum.  The 
results  in  Fig.  3  therefore  provide  ample  evidence  that  ab- 

#  sorption  of  light  by  excited  ioos  does  not  occur  and  that  the 
present  fourfold  laser  operates  by  cooperative  means. 

This  argument  also  rules  out  mechanisms  combining 
cooperative  and  ESA  processes.  For  example,  ESA  of  l.S 
/im  light  fiom  ^S^/2  or  ^Hn/i  states  populated  by  trio  up- 
converskm^^  could  in  principle  reach  the  high  lying 

#  state  leqxHisiUe  for  laser  action  in  this  woric,  but  such 
transitions  are  off  resonant  by  over  218  cm'"^  in  LiYF4  (or 
all  Stalk  levels^^  and  are  not  present  in  the  excitation  spec¬ 
trum.  ESA  transitions  from  ^S^/2  to  ^G^ii/2  are  stiU  farther 
off  resonance.  On  the  other  hand,  quartet  upconversion  can 
im)vide  very  nearly  resonant  excitation  of  the  state 

#  lying  immediately  above  the  upper  laser  level  [generating 

the  fluorescence  at  380  nm  in  Hg.  1(a)].  The  calculated 
energy  defects  on  the  quartet  transitkms 
0,3,3A)^*G  „/2(l)  and  for 

example,  are  only  4  and  7  cm  *  in  LiYF*,  respectively. 
The  bracketed  arguments  ol  the  initial  state  indicate  Stark 

#  levels  of  the  four  interacting  atoms,  and  the  single  argu¬ 
ment  of  the  final  state  gives  the  Stark  levd  of  the  accqitor. 
Taken  together  with  the  close  match  between  the  exdte- 
tion  spectrum  and  ground  state  absorption  (Fig.  3),  this 
argues  strongly  against  any  contributimis  by  ESA  pro- 


FIO.  3.  Infrared  excitation  spectrum  of  fixirfold  upconversion  laser  emis¬ 
sion  at  70I.S  am  ( 7*— 10  K,  lower  trace)  and,  for  comparison,  the  Er^'^ 
abimption  spectrum  (T=9  K,  inverted). 


cesses  to  ^H^/2  population  when  l.S  ftm  pimping  is  used. 

The  implication  that  laser  operation  may  be  sustained 
by  a  cooperative  quartet  process  raises  fundamental  ques¬ 
tions  as  to  just  how  effective  higher  order  cooperative  up¬ 
conversion  processes  can  be  m  general  In  this  regard,  we 
wish  to  point  out  that  an  enhancement  mechanism  exists 
for  cooperative  upconversion  ^ocesses  which  can  signifi¬ 
cantly  extend  their  usefulness  compared  to  multiphoton 
absorption.  For  emission  terminating  on  or  near  levds  in 
which  strmg  cooperative  interactions  occur,  no  effitgney 
“penalty”  b  incurred  for  m-fedd  tqicoaveraiQft  whe^  the 
integer  m  may  be  arbitrarily  large.  When  tqiconvoiHcni 
emission  terminates  on  an  energy  levd  which  can  renew 
cooperative  upconversion,  excited  aUms  are  recycled  with¬ 
out  external  pumping  to  enhance  quantum  eflfciency  by  a 
factor  of  ( 1  —V\ViV3Vmp)~\  where  17^  is  the  m-fohl  co(q)- 
erative  upomversion  efficiency.  i||,  and  1^3  are  Iwanch- 
ing  ratios  for  the  decay  processes  imlicated  in  Hg.  4.  En¬ 
ergy  conversion  between  input  wavdength  and  output 

at  can  be  much  higher  than  the  nruuhniiin  value  n _ 

100/m)%  based  on  a  picture  in  which  only  one 
of  m  atoms  is  upconverted  foDbwing  the  absorptkm  of  m 
pump  photons.  Ignoring  decay  to  dm  ground  state,  the 
energy  efficiency  tdus  on  die  enhanced  value 

(i) 

and  lit  >  i|aMx  when  iqwonverted  atoms  reside  kmg  enough 
in  their  Iwest  excited  stale  to  participate  repeatedly  in 
upconversion.  With  recydin^  energy  efficiency  can  ap¬ 
proach  100%  [for  i||=®l:  (i>sO(:Li4l  widil  a  neghgiUe  en¬ 
ergy  defect  such  that 

In  practice^  dieoreticid  eifonnceni^^qt  ipiantum  effi¬ 
ciency  by  cooperative  m^^  by 

losses  due  to  an  iaetcasiad  ainiiij^  qfl^  (radia¬ 

tive)  levds  on  real  atom  wt^ioBt^^  Branching 
ratios  are  *yr**^^v¥^  tkai^iil|^^pP:d||ill  a  statistical  de¬ 
crease  is  to  be  09^^  “dus¬ 

ters’*  with  appioptia^:h|iii^ai|4#'siltiil^  less  than  a 
nominal  ended  ratfim-  BiMfotar,  ttn  inain  conclusion  is 
that  a  recycling  mechanism  exists  wherdiy  cooperative  dy¬ 
namics  can  mediate  unexpectedly  efficient  upconversion 
emission  for  arbitrarily  high  degrees  (m)  of  upconversion. 
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This  has  important  implications  not  only  for  upcon- 
version  lasers,  but  also  for  other  lasers  in  which  coopera¬ 
tive  nonlinear  dynamics  occur  due  to  heavy  rare-earth  dop¬ 
ing.  One  example  is  the  2.8  fim  Er  laser  reported  by 
Stoneman  and  ^terowitz.'*  These  authors  obtained  an  ef¬ 
ficiency  in  practice  which  was  marginally  greater  than  the 
maximum  theoretical  quantum  efficiency  (ignoring  coop¬ 
erative  dynamics).  A  escalation  of  enhanced  quantum  ef¬ 
ficiency  using  the  formula  above,  ignoring  spontaneous  ra¬ 
diative  losses  for  simplicity,  shows  that  for  980  nm 
pumping,  complete  recycling  of  terminal  level  population 
through  pair  upconversion  results  in  a  “renormalized” 
quantum  efficiency  of  2  (instead  of  1)  and  a  maximum 
energy  efficiency  of  70%  (instead  of  35% ) .  Similar  but  less 
dramatic  improvements  due  to  cooperative  enhancement 
of  quantum  efficiency  are  expect^  for  2.8  fira  lasers 
pumped  by  other  wavelengths.’  Cooperative  enhancement 
probably  also  accounts  for  the  remarkable  efficiency  of  the 
“self-terminating”  Er  trio  laser  at  855  nm." 

In  the  present  experiment,  we  note  that  the  701.5  nm 
upconversion  laser  line  terminates  on  the  *Iu/2  level  of 
trivalent  erbium.  This  level  is  known’’  to  be  effective  in 
pair  upconversion  to  the  *F^/2  level,  but  is  not  the  cooper¬ 
ative  pumping  level  responsible  for  our  inversion.  The 
*Ii3/2  level,  not  the  level,  is  responsible  for  coopera¬ 
tive  pumping  when  excitation  is  tuned  to  1.5  /im.’’"  This 
situation  therefore  resembles,  though  it  is  not  identical  to, 
the  general  enhancement  scheme  of  Fig.  4,  and  we  suggest 
that  enhanced  upconversion  efficiency  through  cooperative 
recycling  may  help  account  for  successful  operation  of  this 
laser,  with  its  unprecedented  (fourfold)  degree  of  upcon¬ 
version. 

In  summary,  we  have  demonstrated  operation  of  an 
upconversion  laser  in  which  population  inversion  is  sus¬ 
tained  in  an  upper  laser  level  at  nearly  four  times  the  in¬ 
cident  photon  energy,”  through  a  fourfold  upconversion 
process.  The  upconversion  mechanism  is  thought  to  be  ex¬ 
clusively  cooperative  in  nature,  arising  from  an  atomic 
“quartet”  interaction  quite  distinct  from  sequential  energy 
transfer"  or  multiphoton  absorption.  The  general  notion  of 
recycling  of  excited  state  population  has  been  introduced 
and  discussed  as  a  possible  mechanism  for  enhancing  mth 


order  upconversion  processes  and  accounting  for  success¬ 
ful  operation  of  the  present  fourfold  upconversion  laser. 
Extension  of  this  result  to  fourfold  laser  operation  at  410 
nm,  555  nm,  and  other  wavelengths  appears  feasible.  In  the 
case  of  the  555  nm  transition,  the  terminal  level  would  be 
in  the  *Ii3/2  manifold,  and  especially  strong  cooperative 
enhancement  of  the  efficiency  may  be  possible. 
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quential  pair  processes  and  sequential  muhiplicrton  absorption  processes 
are  ruled  out  not  only  by  our  previous  worit  on  time  dependence  of 
upconversion  fluorescence,  but  also  by  the  exciution  spectra  which  are 
not  convolutions  of  ground  and  excited  state  resonances. 

*^F.  Auzd.  J.  Lumin.  45,  341  (1990). 

’^R.  C  Stoneman  and  L.  Esterowitz,  Opt  Lett  17,  816  (1992). 

Chou,  Ph.D.  diaaertation,  Massachusetts  Institute  of  Technology, 
Tech.  Rep.  26,  February  1989. 

^^A  pfdimhiary  report  of  these  results  was  flrst  presented  by  S.  C  R.  at 
the  Annual  Meeting  of  the  Optical  Society  of  America,  AHMiquerque, 
New  Mexico,  Sept  20-25  (1992),  mvited  p^per  MJ2.  Results  at  ditfer- 
ent  output  wavelengths  have  been  reported  by  R.  M.  Macfarlane,  E.  A. 
Whittaker,  and  W.  Lcnth,  Eketron.  Lett  28,  2136  (1992). 
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ABSTRACT 

We  report  a  new  mechanism  of  intrinsic  bistability,  operative  in  excited 
atomic  systems  with  substantial  atom-atom  coupling.  Predictions  and 
observations  of  bistable  luminescence  from  Yb3+  pairs  in  crystalline 
Cs3Y2Br9:  Yb3+ are  in  good  agreement  and  have  interesting  implications  for 
local  field  effects  in  solids  and  dense,  ultracold  gases. 
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As  first  reported  by  Gibbs  [1]  many  years  ago,  the  transmission  of 
nonlinear  optical  systems  can  sometimes  exhibit  two  stable  values  for  a 
single  input  intensity,  a  phenomenon  known  as  bistability.  Because  of 
potential  applications  in  all-optical  switching  for  communications  and 
optical  computing,  many  different  mechanisms  of  bistability  have  been 
studied.  Mirrorless  or  intrinsic  bistability  in  dense  media  was  first 
considered  by  Bowden  and  Sung  [2]  and  elaborated  by  subsequent  authors 
[3-6]  who  found  that  ground  state,  near  dipole-dipole  (NDD)  interactions 
responsible  for  making,  the  local  or  Lorentz  field  different  from  the 
incident  field  could  cause  optical  switching.  To  our  knowledge  however 
there  have  been  no  experimental  reports  of  all-optical  switching  medidted 
by  NDD  interactions  nor  any  theoretical  treatments  considering  resonant 
contributions  from  cooperative  excited  state  dynamics.  Here  we  extend 
the  two-level  density  matrix  model  to  incorporate  new  contributions 
expected  from  strong,  dipolar  interactions  among  excited  states  [7] 
pertinent  to  a  variety  of  systems,  including  quasi-resonant  energy 
exchange  collisions  in  gases  [8],  avalanche  upconversion  [9]  and 
cooperative  upconversion  [10]  in  solids.  Intrinsic  bistability  dependent  on 
a  combination  of  ground  and  excited  state  couplings  is  predicted  and  this 
phenomenon  observed  in  cooperative  pair  luminescence  of  the  dimeric 
compound  Cs3Y2Brg:Yb3+  for  the  first  time. 

Pairs  of  atoms  in  dense  systems,  whether  gases  or  solids,  can  relax 
through  energy  transfer  mediated  by  exchange  or  multipolar 
electromagnetic  interactions  at  close  range,  as  depicted  in  Figures  1(a) 
and  1(b).  In  dense  Cesium  vapor  for  example,  an  atom  prepared  in  the  6d 
^D3/2  level  can  relax  non-radiatively  to  the  6p  2P3/2or2Pi/2  level  through 
a  collision  with  a  ground  state  atom  which  undergoes  a  simultaneous, 


energy-conserving  dipole  transition  upward  to  6p  [8].  Such  a  process 
corresponds  to  Fig.  1(a).  As  a  second  example,  in  crystals  doped  with 
Tm3+  impurities,  ions  prepared  in  the  4  state  may  decay  non-radiatively 
to  3F4  if  a  neighbor  makes  a  simultaneous  transition  to  3F4  in  a  condensed 
matter  process  [11]  referred  to  as  cross  relaxation.  In  either  example, 
excitation  energy  is  partitioned  between  the  excited  atom  and  a  ground 
state  "collision"  partner  at  a  rate  which  is  bilinear  in  their  respective 
densities,  and  therefore  nonlinear  with  respect  to  occupation 
probabilities.  Somewhat  surprisingly,  under  selected  initial  conditions, 
nonlinear  dynamics  of  this  kind  can  dominate  collisional  relaxation  in 
gases  [12]  or  lead  to  cross  relaxation  avalanches  in  solids  [11].  The 
reverse  process,  sketched  in  Fig.  1  (b)  and  called  cooperative  upconversion, 
can  also  be  efficient  and  is  the  process  responsible  for  the  observations 
reported  here. 

Our  upconversion  experiments  were  performed  in  single  crystals  of 
Cs3Y2Br9:10%Yb3+,  one  of  a  family  of  compounds  CS3M2X9  (M=Ho3+...Lu3+,  Y3+ 
when  X*CI':  M=Sm3+...Lu3+,  Y3+  when  X-Br)  crystallizing  with  space  group 
R3c  such  that  rare  earth  dimer  units  [M2X9]3-  form  with  trigonal  axes 
coinciding  with  the  trigonal  crystal  axis.  A  key  feature  of  these  dimer  or 
pair-forming  materials  is  that  the  M-M  distance  between  lanthanide  ions 
is  much  shorter  than  in  most  compounds  [13],  including  the  phosphate 
crystal  in  which  Yb3+  cooperative  luminescence  was  observed  originally 
[14].  Since  the  coupling  between  lanthanides  depends  strongly  on  inter¬ 
ion  separation,  materials  of  this  type  exhibit  intense  cooperative 
upconversion  emission,  and  give  other  indications  of  strong  ground  and 
excited  state  interactions.  Neutron  scattering  measurements  of  the  Yb^* 


ground  state  electronic  splitting  in  Cs3Yb2Br9  [15]  for  example  yield  the 
value  3.0  cm-^,  which  is  exceptionally  large  for  rare  earth  ions.  A  low 
power,  continuous-wave  Ti:Al203  laser  tuned  to  various  Stark  components 
of  the  2F7/2  ^^5/2  transition  of  Yb3+  was  adequate  to  excite  green  pair 

luminescence  at  a  wavelength  of  500  nm  in  samples  cooled  to  liquid 
helium  temperatures.  The  intensity  of  the  emission  was  carefully 
recorded  as  temperature  and  incident  intensity  were  slowly  varied. 

The  dynamics  of  interest  (Fig.  1(b))  can  be  treated  with  a  two-level 
pair  model  of  cooperative  processes  (Fig.  2).  This  model  incorporates 
ground  state  dipolar  interactions  which  produce  local  field  corrections  in 
the  optical  polarization,  following  previous  work  [5,  16].  However, 
excited  state  interactions  are  also  included.  Of.  the  three  states  shown  in 
Fig.  1,  we  note  that  one  participates  merely  as  an  intermediate  state.  Its 
role  as  a  virtual  state  can  nevertheless  be  preserved  in  a  two-level  model, 
as  illustrated  in  Fig.  2(a)  and  2(b),  by  the  incorporation  of  relaxation 
terms  in  the  equations  of  motion  with  density-dependent  rates.  For 
example  in  Fig.  2(b),  cooperative  upconversion  mediates  decay  of  state  |2> 
population  P22  at  a  rate  of  2a.p22^,  where  a  is  the  pair  upconversion  rate 
constant  and  the  factor  of  2  refers  to  the  two  excited  atoms  which  decay 
for  each  cooperative  event.  This  contributes  a  term  -<x(ff3-ff3o)2  to  decay 
of  the  inversion  ff3*P22”Pii  with  respect  to  its  equilibrium  value  cr3o  in  the 
absence  of  fields,  provided  there  is  negligible  thermal  occupation  of  the 
excited  state  (ffso—l)-  From  the  Schrodinger  equation  of  motion  for  the 
density  matrix,  the  dynamic  equations  for  the  off-diagonal  coherence  a  21 
and  the  inversion  <73  are  therefore: 
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(1) 


=  -i  (A  +  eoj)  a^i  - 1  iQc^  -  Ta^, 

03  =  -i  a  (Gjj  -  0*21)  -  7(a3  -  G30)  -  0(03  -  03/  (2) 


Here  c=n|j,2/(sfico)  is  the  local  field  correction  term,  which  has  units  of 
frequency,  just  like  a.  n  is  the  number  density,  is  the  dipole  transition 
moment,  s  is  a  factor  (equal  to  3  in  cubic  solids)  determined  by  a  lattice 
sum  [17],  CO  is  the  permittivity  of  vacuum,  y  is  the  excited  state  decay 
rate  and  r  the  dephasing  rate.  The  detuning  of  light  from  the  resonance 
frequency  at  uq  is  a=«o-w,  and  is  the  Rabi  frequency.  cr2i=-ip2i  is  a 

convenient,  slowiv-varying  amplitude  of  the  off-diagonal  density  matrix 
element  p2i*P2iexp[-i(Qt-k.r)].  Overdots  indicate  time  derivatives.  ' 
From  (1)  and  (2)  the  population  inversion  <13  and  the  slowly- varying 
envelope  of  the  optical  polarization  P=  2nipi2cr2i  niay  be  calculated.  Since 
cooperative  upconversion  and  cross  relaxation  are  inverse  processes,  they 
are  both  described  by  very  similar  dynamical  equations  and  conclusions 
drawn  below  also  pertain  to  avalanche  upconversion  processes  (Figs.  1(a), 
2(a)).  In  steady-state  conditions  the  time  derivatives  in  (1)  and  (2)  equal 
zero.  This  yields 


-yiQa3[r-i(A  +  eo3)] 

a2i=-"^— ^ - 2 - 

r +  (A  +  eo3)^ 

By  substituting  (3)  into  (2)  and  simplifying,  we  find 

ao^  +b<^  +  co^  +  dO3  +  e  =  0 

where  the  coefficients  are  given  by 


(3) 


(4) 
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a  =  ae^,  b  =  c  =  -Y03Qe^+2YAe+ao^Qe^+ar^-4aa3QAe+aA^, 

d  =  -2YCT3QAe+Q^r+Yr^ +2ao^QAe-2ao3Qr^+YA^-2aa3QA^,  and 
e  =  ao2/^-Y03or^-Y03oA^  +  ao^^A^  . 


The  quartic  equation  in  (4)  has  solutions  which  are  multi-valued  and 
describe  the  full  features  of  bistability  to  be  expected  in  the  inversion  of 
cooperative  dynamic  systems.  When  the  upconversion  rate  a  is  zero,  the 
coefficient  of  the  quartic  term  cxc2  vanishes  and  we  recover  the  results  of 
Refs.  [2-6,16].  However  when  both  cx  and  c  are  not  zero,  the  fourth  order 
coefficient  becomes  finite  and  leads  to  changes  in  the  nonlinear  dynamics 
which  are  the  joint  result  of  ground  and  excited  state  NDD  interactions.* 

In  Fig.  3,  the  solution  of  (4)  for  the  steady-  state  inversion  as  is  used 
to  plot  (ff3-a3o)^.  a  quantity  proportional  to  the  pair  luminescence 
intensity,  as  a  function  of  incident  intensity  Without  upconversion, 
the  leftmost  curve  shows  that  weak  bistability  is  predicted  when  the 
local  field  factor  c  exceeds  excited  state  relaxation  frequencies  y  and  r. 
The  solutions  in  the  central  portion  of  this  plot  between  the  turning 
points  (dashed)  are  demonstrably  unstable,  causing  hysteresis  in  the 
inversion  versus  incident  intensity.  As  intensity  is  increased  from  zero, 
the  predicted  pair  luminescence  intensity  at  first  follows  the  lower 
branch  of  the  curve,  but  reaches  a  point  at  which  the  only  stable  solution 
available  at  higher  powers  lies  on  the  upper  branch  of  the  curve.  Here, 
-switching  to  the  upper  branch  occurs  and  switch-down  only  occurs 
subsequently  if  the  incident  intensity  falls  below  the  point  at  which  only 
a  lower  branch  solution  exists. 

As  the  upconversion  rate  a  is  increased  from  zero  to  a  value 
comparable  to  c  (curve  farthest  to  the  right  in  Fig.  3),  two  significant 


effects  occur.  First,  a  higher  input  power  is  required  to  reach  bistability. 
Second,  an  immense  magnification  of  the  area  of  the  hysteresis  loop 
occurs  with  only  minor  changes  in  the  relative  shape  of  the  curve.  This  is 
evidence  that  excited  state  dipole-dipole  interactions  influence  nonlinear 
dynamics  in  much  the  same  way  ground  state  interactions  do,  in  keeping 
with  the  fact  that  both  enter  the  coefficient  of  the  quartic  term  on  a 
similar  footing  as  the  product  ac2.  However,  for  reasons  not  yet 
completely  understood  but  possibly  related  to  the  resonant  nature  of 
excited  state  interactions  considered  here,  the  effect  of  a  is  much  greater 
than  that  of  c. 

In  Fig.  4,  experimental  results  are  shown  which  confirm  basic 
predictions  of  the  model.  In  Fig.  4(a)  upconversion  emission  abruptly 
increases  at  a  well-defined  intensity  as  the  incident  power  increases. 
Sample  emission  maintains  an  almost  constant  value  as  the  incident 
intensity  is  increased  further  and  switches  down  only  when  the  power  is 
decreased  well  below  the  original  switching  point.  As  the  atom-atom 
coupling  increases  toward  lower  temperatures,  the  hysteresis  loops 
enlarge  and  move  to  higher  excitation  energies  as  expected.  However  the 
step  height  enlarges,  in  a  trend  which  is  not  understood  at  present.  This 
may  result  from  neglect  of  resonant  energy  migration  in  the  excited  state, 
which  cannot  be  included  using  single  atom  basis  states  in  a  2-level 
model  although  It  would  undoubtedly  enhance  the  importance  of 
cooperative  interactions.  In  Fig,  4(b),  similar  hysteresis  is  apparent  when 
the  incident  intensity  Is  held  fixed  and  temperature  is  varied.  In  this  case 
the  effective  inter-atom  coupling  changes,  because  the  overlap  integral 
and  competing  decay  processes  are  temperature  dependent.  Hence  this  is 
roughly  equivalent  to  varying  a  at  a  fixed  intensity  in  our  model,  in  which 


case  the  theory  predicts  very  similar  bistability  results.  Both  sets  of 
observations  clearly  show  that  cooperative  emission  exhibits  hysteresis 
without  feedback  (mirrors). 

Broad  inferences  may  be  drawn  from  the  present  work  by  adopting 
the  following  perspective.  In  both  gases  and  solids,  cooperative 
interactions  like  cross  relaxation  or  upconversion  may  be  pictured  as 
"collisional"  processes  which  occur  on  timescales  comparable  to  or 
shorter  than  the  average  collision  time.  In  gases,  collision  times  are 
typically  very  short,  and  collisional  interactions  brief.  On  the  other  hand, 
near-neighbor  rare  earth  impurities  in  solids,  separated  by  short 
distances  as  if  in  a  "frozen"  collision,  interact  in  excited  states  on 
timescales  limited  by  lifetimes  of  the  energy  levels,  which  can  be  as  long 
as  ten  milliseconds.  This  analogy  with  gaseous  collisions  is  in  fact  more 
apt  than  generally  recognized,  since  collisional  depolarization  ratios  of 
rare  gases  for  example  can  be  calculated  from  elasto-optic  coefficients 
of  rare  gas  solids  [18].  Hence  when  long  interaction  times  are 
encountered,  such  as  in  rare  earth  solids  or  for  nearly  stationary  atoms  in 
optical  traps,  the  concept  of  slow  collisions  can  be  invoked  to  unify  the 
seemingly  distinct  domains  of  coupled  dopants  in  solids  and  colliding 
ultracold  atoms  in  vacuum.  The  model  presented  here,  in  which  atoms  are 
assumed  to  be  effectively  at  rest,  is  expected  to  encompass  both  physical 
situations. 

Although  not  explicitly  reported  here,  optical  absorption  in  this  and 
other  coupled  systems  is  also  expected  to  undergo  switching,  since  it  is 
proportional  to  53.  That  is,  absorption  by  coupled  atoms  In  dense  excited 
systems  with  transition  energy  coincidences  should  exhibit  hysteresis. 
This  implies  that  many  related  effects  may  be  observable  in  dense  vapors 


and  gases.  For  example,  the  absorption  behavior  of  ultracold  colliding 
atoms  which  happen  to  have  nearly  identical  ground  and  excited  state 
transition  energies  should  exhibit  hysteresis  on  selected  transitions.  At 
microkelvin  temperatures  and  below,  atoms  are  essentially  at  rest  on  the 
time  scale  of  excited  state  relaxation  processes.  Consequently,  we  expect 
behavior  related  to  that  calculated  and  observed  here  for  stationary  atoms 
in  a  concentrated  rare  earth  solid  to  apply  to  ultracold  dense  gases  too. 

In  summary,  we  have  shown  that  pair  emission  of  Yb^^  ions  in 
Cs3Y2Br9:Yb3+ exhibits  intrinsic  bistability  as  a  function  of  incident 
intensity  and  temperature.  A  model  in  which  contributions  to  the 
nonlinear  polarization  from  ground  and  excited  state  interactions  between 
Yb3+  ions  are  treated  on  an  equal  footing  gives  good  agreement  with 
experiments.  Furthermore,  the  mere  fact  that  the  first  example  of 
bistability  due  to  atom-atom  interactions  has  emerged  from  a  cooperative 
upconversion  system  supports  the  general  conclusion  that  resonant 
excited  state  interactions  amplify  nonlinear  response  associated  with  the 
Lorentz  field,  rendering  bistability  much  easier  to  observe.  It  confirms 
experimentally  that  cooperative  nonlinearities  can  induce  hysteretic 
behavior  in  crystals  without  a  cavity.  Bistability  from  this  mechanism 
can  be  expected  to  occur  also  in  avalanche  systems  which  decay  by  run¬ 
away  cross  relaxation.  Furthermore,  while  our  results  are  most  pertinent 
to  stationary  atoms  in  solids,  we  have  argued  that  they  have  broader 
implications  for  selected  transitions  in  vapors  which  induce  quasi¬ 
resonant  collisional  interactions.  Similar  effects  may  occur  in  high 
density  optical  traps  below  the  Doppler  cooling  limit  when  internuclear 
separation  does  not  change  appreciably  during  an  excited  state  lifetime. 


Finally,  these  results  suggest  new  extensions  of  theories  of  the  local 
field  [17]  in  dense  atomic  systems. 
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FIGURE  CAPTIONS 


Figure  1.  Cooperative  dynamics  involving  two  3-level  atoms,  giving  rise 
to  (a)  cross  relaxation,  and  (b)  pair  upconversion.  Wiggly  arrows  indicate 
photons.  Curved  arrows  indicate  (non-radiative)  coupled  atom 
(cooperative)  dynamics. 

Figure  2.  Two-level  theoretical  models  of  the  cooperative  dynamics 
depicted  in  Fig.  1  in  which  the  third  (intermediate)  states  are  treated  as 
virtual  states.  Wiggly  arrows  indicate  photons.  Curved  arrows  indicate 
(non-radiative)  cooperative  relaxation  processes. 

Figure  3.  Calculated  pair  luminescence  intensity  (proportional  to 
(ff3-ff3o)^)  in  a  system  undergoing  pair  emission,  plotted  as  a  function  of 
incident  intensity  (q2)  on  resonance  (A=0).  The  solid  (dashed)  curves 
indicate  stable  (unstable)  solutions.  From  left  to  right,  the  five  plotted 
curves  correspond  to  values  of  inter-atomic  coupling  of  <x=0,  2,  4,  6,  8 
respectively.  All  frequencies  are  referenced  to  y=r=1  and  the  local  field 
parameter  has  the  fixed  value  c=10  ccr,  where  ccr  »3V3  is  a  critical  field 
needed  for  bistability  when  a=0. 

Figure  4.  Upconversion  emission  intensity  at  A=  500  nm  versus  incident 
laser  power  from  Yb3+-Yb3+  pairs  in  CssYaBrgilO®/©  Yb3+at  various  fixed 
temperatures  (top  to  bottom):  31 K,  27K,  23K,  19K,  15K,  11K.  (b) 
Upconversion  emission  intensity  at  500  nm  versus  temperature  in 
Cs3Y2Br9:10%  Yb3+ at  various  fixed  incident  intensities  (top  to  bottom): 
814  ,  640.  539,  415  W/cm2. 
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A  INTRODUCTION 

Many  of  the  properties  which  make  diamond  uniquely  attractive  for  mechanical  and  electronic 
applications  also  lend  interest  to  its  potential  for  optical  applications.  In  particular,  its  wide 
bandgap  and  the  exceptional  stability  of  various  optical  centres  in  this  host  make  it  attractive 
for  generation  of  radiation  throughout  the  visible  spectral  region.  Furthermore,  the  high 
thermal  conductivity,  excellent  hardness,  and  high  damage  threshold  of  diamond  suit  it  weU  to 
high  power  laser  operation.  Here  we  present  a  review  of  progress  and  prospects  in  diamond 
laser  development 


In  Section  B,  a  few  radiation-induced  defect  centres  which  may  be  useful  for  the  generation  of 
laser  radiation  in  diamond  are  compared  with  high  gain  laser  centres  of  the  alkali  halides.  In 
Sections  C  and  D,  progress  on  H3  laser  development  and  prospects  for  other  diamond  centres 
are  discussed. 


B  OPTICAL  CENTRES  FOR  LASERS  IN  DIAMOND 

Pure  diamond  is  an  indirect  gap  material  like  silicon.  As  a  consequence,  there  is  little  hope  of 
obtaining  stimulated  emission  directly  at  the  bandgap  recombination  wavelength  of  227  nm 
[1].  However  intrinsic  point  defects,  impurities,  and  centres  composed  of  both  point  defects 
and  impurities  provide  potential  candidates  for  laser  action  at  longer  wavelengths.  Here,  we 
tabulate  spectroscopic  parameters  and  theoretical  laser  gain  for  a  few  centres  which  combine 
impurities  and  vacancies  and  review  their  potential  for  laser  applications.  Procedures  for  the 
synthesis  of  centres  and  general  characteristics  of  colour  centre  lasers  may  be  found  in  [2]  and 
[3]  respectively. 

From  TABLE  1  two  general  conclusions  may  be  drawn.  First,  diamond  colour  centres  have 
exceptionally  broad  fluorescence  linewidths  Av  compared  to  those  of  the  commercially 
successful  alkali  halide  lasers,  exemplified  by  the  F,"'  laser.  As  a  result,  their  potential  tuning 
range  is  many  times  greater.  Second,  their  gain  coefficients  are  lower  than  that  of  Fj^'rKCl, 
but  still  very  respectable  considering  the  breadth  of  emission.  The  order  of  magnitude 
^ffetence  between  F,"^  and  N-V  gain  coefficients  may  be  overcome  in  principle  by  merely 
increasing  the  excited  state  density  by  a  factor  of  ten. 
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TABLE  1  Comparison  of  diamond  and  F2*;KC1  colour  centre  gain  coefficients 
yfor  constant  inversion  density  N  =  lO'*  cm'^. 


n 

T 

(K) 

"^11 

(ns) 

n 

1.4841  [3] 

77 

Tj/n  =  80  [3] 
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2.4065  [4] 

295 

13.3  [5] 

0.99  [5,6] 

2.4262  [4] 

295 

16  [2] 

0.95  [7] 

(10‘^  s') 


1.69  [3 
6.53  [5 
5.27  [5 


I 


^0 

(pm) 

Y 
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1.680  [3] 

3.52 

0.697  [5] 

0.366 

0.531  [5] 

0.201 

Since  the  discovery  of  the  first  diamond  laser  [8],  two  other  experiments  have  demonstrated 
laser  action  in  diamond  using  H3  colour  centres.  In  this  unpublished  work,  high  densities  of 
H3  centres  in  natural  and  synthetic  diamonds  were  used,  as  described  briefly  below. 

In  1988  at  the  Itami  Research  Laboratories  of  Sumitomo  Electric  in  Japan,  a 
flashlamp-pumped  Coumarin  dye  laser  with  pulsewidth  of  500  ns  was  used  to  side-pump  a 
natural  diamond  crystal  of  dimensions  3x2x2  mm^  at  an  excitation  wavelength  of  490  nm 
inside  a  2-mirror  resonant  cavity.  The  end  mirror  was  a  total  reflector  with  a  radius  of  40  cm 
and  the  output  mirror  had  the  same  curvature  and  a  reflectivity  of  97%.  Both  mirrors  were 
spherical  and  separated  by  approximately  35  cm.  Laser  action  in  the  yellow  spectral  region 
was  obtained  at  room  temperature  with  a  threshold  of  3  MW/cm"  without  auxiliary  cooling 
[9].  Laser  output  ceased  when  cavity  mirrors  were  detuned  or  blocked,  thus  extending  the 
original  result  by  the  use  of  an  external  cavity. 

Continuous- wave  laser  action  was  demonstrated  in  the  same  year  at  the  University  of 
Strathclyde  [10].  Optical  excitation  was  performed  with  a  continuous-wave  argon  ion  laser 
operating  at  488  nm.  In  this  case  a  3-mirror  cavity  was  used  to  compensate  for  astigmatism 
introduced  by  tilting  the  sample  to  Brewster’s  angle  for  minimum  reflective  losses,  and  the 
crystal  was  cooled  to  77  K.  Output  power  is  shown  as  a  function  of  input  power  in  FIGURE 
1  for  two  synthetic  diamond  samples  from  De  Beers.  Although  slope  efficiency  was  very  low 
(6.6  X  10"^%),  the  tuning  range  of  this  laser  was  remarkably  wide,  comparable  with  the 
fluorescence  spectrum  shown  in 
FIGURE  2.  Its  low  efficiency 
may  be  related  to  residual  N-V  io 

absorption  in  the  laser  emission 
region,  or  to  transient  losses  ^ 

associated  with  B  or  Ni  3 

impurities  used  to  facilitate  H3  “S  5 
formation  during  crystal  growth.  3 
Sumitomo  researchers  have  also  ® 
synthesised  H3  crystals  [11]  with 
very  bright  photoluminescence,  q 
but  to  date  have  not  reported  0  1.0  2.0 

continuous  lasing.  Input  Powor  (W) 

FIGURE  I  Coniinuous-wave  laser  output  power  as  a  function  oi 

input  pump  power  at  488  nm  in  two  synthetic  diamond  crystals 

(De  Beers  JH5.  RDl)  at  T  =  77  K.  (With  permission  [10].) 
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FIGURE  2  Comparison  between  tuning  range  of  the  CW  diamond  laser  (solid  circles) 
and  the  fluorescence  spectrum  (solid  curve)  of  the  laser  active  H3  centre,  obtained 
by  time  resolved  photoluminescence.  (With  permission  [10].) 

D  PROSPECTS  FOR  NEW  DIAMOND  LASERS 

The  H3  laser  described  above  is  the  only  laser  of  any  kind  able  to  scan  continuously  through 
the  wide  visible  range  from  500  -  600  nm.  This  unique  capability  argues  strongly  for  renewed 
research  to  improve  crystal  synthesis  in  an  effort  to  minimise  parasitic  losses.  Additionally 
however,  other  promising  approaches  for  lasers  at  both  shorter  and  longer  wavelengths  in 
diamond  should  be  considered.  Several  of  these  are  outlined  below. 

The  theoretical  gain  coefficient  for  the  N-V  centre  (TABLE  1)  is  nearly  twice  that  of  the  H3 
centre.  Hence  it  is  an  excellent  laser  candidate,  compatible  with  efficient  diode  laser  pumping 
schemes.  Of  course,  absoiptive  losses  in  the  emission  band  must  be  kept  very  low  for  net 
positive  gain  to  be  realised  experimentally.  Encouragingly,  the  only  centre  with  broadband 
absorption  in  the  range  650- 750  nm  is  the  GRl  vacancy  centre,  which  exhibits  ver>  ow 
absorption  (and  emission)  at  room  temperature  and  can  be  removed  by  appropriate  annealing 
procedures.  Also,  the  intrinsic  quantum  efficiency  of  the  N-V  centre  is  unity  [6].  Although 
recent  spectroscopic  studies  [12]  confirmed  the  presence  of  a  metastable  level  below  the 
638  nm  fluorescent  level  of  the  N-V  centre  which  acts  as  a  sink  for  excited  state  population, 
its  effect  is  negligible.  Furthermore,  high  concentrations  of  this  centre  can  be  produced  by  a 
simple  procedure,  so  its  overall  prospects  for  laser  action  remain  good. 

Other  nitrogen-related  centres  may  be  viable  laser  candidates.  Scant  information  is  available 
on  detailed  models  or  stability  of  such  centres,  however,  and  few  have  strong  broadband 
emission  in  the  visible  or  near  infrared  (IR)  spectral  regions.  Other  approaches  nevertheless 
offer  promising  new  directions.  For  example,  changing  diamond  into  a  direct  gap  material 
through  the  use  of  strained  layer  structures  has  been  considered  [13].  Additionally, 
radiation-induced  and  aggregate  impurity  centres  in  boron-doped  diamond  should  be 
investigated.  Boron  is  easily  substituted  into  the  diamond  matrix  and  results  in  useful 
semiconductivity.  By  incorporating  other  impurities  which  fit  reasonably  well  into  the  lattice 
(C  ionic  radius  =  0.16  A),  like  nitrogen,  sulphur,  chlorine,  phosphorus  and  silicon  (ionic  radii 
of  0.16  A,  0.44  A,  0.37  A.  0.34  A  and  0.42  A  respectively),  many  new  stable  aggregate 
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centres  can  undoubtedly  be  created  and  pumped  electrically,  warranting  further  investigation. 
By  combining  such  centres  with  radiation-induced  vacancies  using  suitable  annealing 
procedures,  additional  possibilities  for  stable  centres  may  be  realised. 

Rate  earth  (RE)  dopants,  which  are  currently  under  investigation  for  lasers  and  display 
components  in  compound  semiconductors  [14],  and  constitute  the  laser  active  species  in 
dozens  of  commercial  solid  state  lasers,  may  also  be  useful  for  semiconductor  diamond  defect 
lasers.  However,  the  ionic  radii  of  trivalent  rare  earth  ions  all  exceed  0.86  A.  Hence  it  seems 
unlikely  that  high  enough  concentrations  of  such  impurities  can  be  achieved  in  diamond  to 
permit  such  applications.  Transition  metal  ions  tend  to  have  smaller  radii  and  may  therefore  be 
more  promising  candidates  for  laser  research  in  this  host  Boron  itself  absorbs  strongly 
throughout  the  visible  spectrum,  so  it  is  likely  that  in  B-doped  semiconducting  diamond  useful 
emission  wavelengths  of  impurity  metal  laser  ions  (possibly  associated  with  compensating, 
radiation-induced  defects)  will  be  restricted  to  the  infrared.  At  very  long  wavelengths  because 
of  its  renowned  transparency  in  the  far  infrared  (FIR),  diamond  might  ultimately  prove, 
however,  to  be  the  material  of  choice,  indeed  periiaps  the  only  material  compatible  with 
solid-state  FIR  laser  operation. 
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Abstract 

Ultrafast  photon  echo  experiments  reveal  highly  modulated  quantum  beats  in  nitrogen-vacancy  centers  in  diamond, 
providing  direct  measurements  of  excited  state  splittings  for  the  first  time.  We  also  confirm  the  existence  of  new  peaks  in 
echo  excitation  spectra  and  report  novel  polarization  dependent  interferences. 


The  establishment  of  definitive  models  of 
luminescent  centers  in  solids,  particularly  in 
widegap  semiconductors,  is  becoming  increasingly 
important  in  the  context  of  next  generation  flat 
panel  displays.  Here,  we  report  the  direct  deter¬ 
mination  of  excited  state  structure  and  dephasing 
dynamics  of  the  nitrogen-vacancy  (N-V)  center  in 
diamond  from  ultrafast  photon  echo  observations. 
These  results  provide  the  basis  for  a  detailed  model 
of  the  center. 

Previous  results  from  isochronal  annealing  ex¬ 
periments  and  uniaxial  stress  measurements  by 
Davies  and  Hamer  [1]  on  the  638  nm  zero  phonon 
transition  of  the  N-V  center  provided  evidence  for 
a  trigonal  center  of  C^v  symmetry  consisting  of 
a  carbon  vacancy  adjacent  to  a  nitrogen  impurity. 
The  ground  and  excited  states  were  determined  to 
be  of  A  and  E  symmetry,  respectively  [1].  Electron 
paramagnetic  resonance  (EPR)  measurements  by 
Loubser  and  Van  Wyk  [2]  confirmed  the  basic 
physical  model  of  the  center.  Transient  hole- 
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burning  experiments  [3]  and  Curie-law  depend¬ 
ence  of  the  magnetic  susceptibility  from  ground 
state  EPR  measurements  at  2.88  GHz  by  Redman 
etal.  [4]  established  the  triplet  character  of  the 
ground  state.  Persistent  hole-burning  measure¬ 
ments  also  revealed  transitions  at  frequencies  com¬ 
bining  the  ground  state  spin-spin  splitting  with 
other  frequencies  on  the  GHz  scale,  providing  in¬ 
direct  evidence  for  ^E  excited  state  structure  [5]. 
Optically  detected  microwave  resonance  (ODMR) 
measurements  [6],  in  which  the  laser  probe  was 
tuned  across  the  ZPL,  also  revealed  a  larger 
40cm” ^  splitting  in  the  ODMR  spectrum,  sugges¬ 
ting  that  an  even  larger  splitting  existed  in  the 
excited  state. 

In  the  present  work,  ultrafast  stimulated  photon 
echo  experiments  were  performed  to  investigate 
excited  state  structure  directly.  A  dual  jet,  syn- 
chronously-pump  DCM  dye  laser  was  employed, 
generating  cavity-dumped  subpicosecond  pulses 
at  a  variable  repetition  rate  near  100  KHz.  The 
sample  was  a  high-quality  synthetic,  nitrogen  rich 
single  crystal.  This  was  irradiated  with  1  7  MeV 
electrons  and  subsequently  annealed  to  produce 


0022-23 13/94/$07.(X)  ©  1994  -  Elsevier  Science  B.V.  All  rights  reserved 
S5£)/ 0022-23  13  (93)  E0263-W 


740 


S.C.  Rand  et  aL  ,  Journal  of  Luminescence  60<&6I  (1994  )  739-  741 


DELAY  (ns) 


Fig.  1.  (a)  Stimulated  photon  echo  decay  data  from  diamond 
N-V  centers,  showing  highly  modulated  quantum  beats  at  6  K, 
X  -  638.2  nm.  Incident  beams  are  all  linearly  polarized  in  the 
<00  1>  direction,  (b)  Power  spectrum  of  echo  data. 


high  concentrations  of  N-V  centers.  Three  input 
beams  were  arranged  in  a  phase  conjugate  geo¬ 
metry  with  counter-propagating  pump  waves  and 
a  probe  wave  incident  at  a  small  angle  with  respect 
to  the  forward  pump.  Pump  waves  were  aligned 
parallel  to  the  [100]  direction.  Spatially  filtered, 
time-averaged  echo  signals  were  recorded  using 
lock-in  detection  techniques  with  a  photomultip¬ 
lier.  Signals  were  recorded  as  a  function  of  probe 
delay,  incident  wavelength,  and  temperature,  and 
they  verified  to  be  independent  of  repetition  rate. 
For  polarization-dependent  studies,  quarter-wave 
plates  were  inserted  into  all  three  incident  beams, 
using  an  additional  quarter-wave  plate  and  polariz¬ 
er  combination  to  analyze  the  resulting  echo. 


The  photon  echo  intensity  is  shown  in  Fig.  1(a) 
as  a  function  of  probe  delay,  revealing  deep  modu¬ 
lations  in  the  photon  echo  decay.  The  wavelength 
was  tuned  to  the  red  side  of  the  ZPL  (/.  =  638.2  nm) 
with  all  incident  beams  polarized  in  the  <001) 
direction.  Figure  1(b)  shows  the  power  spectrum  of 
these  data.  The  exponential  decay  in  Fig.  1(a)  was 
first  divided  from  the  data  before  Fourier  trans¬ 
formation.  We  measured  frequency  components  at 
A  =  2.87  ±  0.01  GHz,  corresponding  to  the  known 
ground  state  splitting,  as  well  as  components  at 
B  =  0.61  ±  0.02  GHz,  C  =  1.65  ±  0.01  GHz,  and 
D  =  2.30  ±  0.01  GHz,  related  to  excited  state  split¬ 
tings.  These  observed  frequencies  agree  with  our 
persistent  hole-burning  observations  [5],  but  in 
addition  confirm  that  all  modulations  originate 
directly  from  splittings  in  the  and  states 
coupled  by  the  light,  rather  than  relaxed  states,  since 
the  signal  appears  as  a  quantum  beat  from  coher¬ 
ence  established  within  the  incident  pulse  duration. 

The  existence  of  two  excited  states  with  differ¬ 
ent  dephasing  behaviors  was  revealed  directly  by 
recording  the  excitation  spectrum  of  the  echo  sig¬ 
nals  across  the  ZPL.  The  excitation  spectrum 
shows  two  distinct  peaks  separated  by  a  width  of 
42±2cm"^  not  resolvable  in  linear  absorption 
measurements  of  these  strongly  inhomogeneously 
broadened  centers.  The  dephasing  rates  versus  tem¬ 
perature  of  these  two  states  were  measured  with  the 
laser  tuned  to  the  red  (638.9  nm)  and  blue 
(636.5  nm)  sides  of  the  ZPL  linear  absorption  peak. 
Dephasing  times  at  5.5  K  for  the  red  and  blue  peaks 
were  4.92  ±  0.24  ns  and  56  ±  8  ps,  respectively,  cor¬ 
responding  to  four  times  the  measured  decay  times 
in  the  photon  echo  signals.  Theoretical  fits  confirm 
that  the  upper  of  the  two  ^E  states  decays  rapidly 
by  a  direct  process  to  the  lower  level  as  one  might 
expect,  while  the  lower  state  dephases  more  slowly 
through  a  virtual  2-phonon  process.  The  dominant 
activation  energy  for  the  direct  process  was 
67  ±  25cm~\  in  reasonable  correspondence  with 
the  measured  splitting  of  the  two  ^E  states.  For  the 
lower  level,  an  Orbach  process  dominates  at  low 
temperatures  with  an  activation  energy  of  39  ± 
10  cm  ■  ^  also  in  good  agreement  with  the  measured 
excited  state  splitting. 

Finally,  a  detailed  investigation  of  the  polariza¬ 
tion  dependence  of  the  echo  intensities  with  circu- 
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larly  polarized  beams  revealed  an  interesting  selec¬ 
tion  rule  which  caused  signal  intensities  to  drop  by 
an  order  of  magnitude  when  the  pulse  polarization 
sequence  was  or  This  effect 

arises  from  interferences  in  the  echo  signal  due  to 
a  superposition  of  third-order  interactions  in  the 
presence  of  ground  state  magnetic  degeneracies. 
For  fixed  probe  delay,  ignoring  quantum  beat  con¬ 
tributions  and  population  relaxation,  the  ith  com¬ 
ponent  of  the  echo  polarization  has  the  form 

Pp'  = 

L®  >  Jkl 


X 


-h 


E^jEu 


(1) 


where  Er,j  is  the  jih  component  of  the  applied  elec¬ 
tric  field  of  the  rjih  pulse,  and  is  a  third- 

order  susceptibility  that  includes  the  tensoral  na¬ 
ture  of  the  interaction.  The  summation  over  a  in 
Eq.  (1)  is  over  the  different  defect  orientations  in  the 
lattice. 

If  we  consider  the  specific  example  ofaa'^a'cr'^ 
sequence  for  a  specific  orientation,  the  x  compon¬ 
ent  of  the  third-order  polarization  in  Eq.  (1)  is 


P.X  =  eo£o|£ol^ 
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(and  equivalent  sequences  with  reversed  circular 
polarization),  respectively.  These  results,  which  will 
be  presented  in  detail  in  a  later  work,  assume  ap¬ 
plied  fields  propagate  in  the  [001]  direction  and 
are  in  good  agreement  with  our  observations. 

In  summary,  we  have  presented  quantum  beat 
measurements  of  excited  state  splittings  in  the  N-V 
center  in  diamond.  The  existence  of  a  large  splitting 
in  the  excited  state  manifold  was  revealed  in  echo 
excitation  spectra.  Temperature-dependent  de¬ 
phasing  studies  of  these  two  states  indicated  that 
direct  and  virtual  2-phonon  processes  were  the 
primary  mechanisms  responsible  for  dephasing 
from  the  upper  and  lower  ^E  levels,  respectively. 
Finally,  we  observed  a  polarization  selection  rule 
which  can  be  explained  as  an  interference  effect  due 
to  Zeeman  coherence  within  the  center  which  is 
consistent  with  our  model  of  the  electronic  struc¬ 
ture. 
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where  Eq  is  the  magnitude  of  the  applied  fields.  In 
a  center  with  cubic  symmetry,  the  only  nonzero 
terms  in  Eq.  (2)  are  xi%y  =  =  xScyW  = 

in  which  case  =  0.  For  the  more  complicated 
case  of  the  C3V  symmetry  N-V  center,  P^  is  still 
reduced  significantly  for  the  case  when  all 

defect  orientations  are  considered.  Explicit  group- 
theoretical  calculations  of  for  a  transition  be¬ 
tween  a  doubly  degenerate  E  symmetry  state  in  the 
ground  state  manifold  and  a  nondegenerate  excited 
state  in  the  N-V  center  predict  relative  echo  inten¬ 
sities  with  ratios  of  64:25:4  for 
or  (T‘^cr”(r“,  and  polarization  sequences 


References 


[1]  G.  Davies  and  M.F.  Hamer,  Proc.  Roy.  Soc.  London  A  348 

■  (1970)  285. 

[2]  J.H.N.  Loubser  and  J.A.  Van  Wyk,  Diamond  Research 
(Industrial  Diamond  Information  Bureau,  London,  1977) 

p.  11, 

[3]  N.R.S.  Reddy,  N.B.  Manson  and  E.R.  Krausz,  i.  Lumin,  38 
(1987)  46, 

[4]  D,A,  Redman,  S.W.  Brown,  R.H,  Sands  and  S.C.  Rand. 
Phys,  Rev.  Lett.  67  (1991)  3420. 

[5]  D.A.  Redman,  S.W.  Brown  and  S.C.  Rand,  J.  Opt.  Soc  Am. 
B  9  (1992)  768. 

[6]  E.  van  Oort,  B.  van  der  Kamp,  R.  Sitters  and  M.  Glasbeek. 
J.  Lumin.  48&49  (1991)  803. 


PHYSICAL  REVIEW  A 


VOLUME  49,  NUMBER  1 


JANUARY  1994 


ARTICLES 


Tunneling  dynamics  of  squeezed  states  in  a  potential  well 
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WKB  calculations  have  revealed  that  tunneling  dynamics  of  localized  wave  packets  in  simple  potential 
wells  depend  sensitively  on  the  degree  of  squeezing  of  the  wave  function.  Squeezed  states  are  shown  to 
tunnel  more  slowly  in  general  than  eigenstates  of  the  well  and  states  which  are  initially  localized  near  the 
bottom  of  a  well  have  much  lower  tunneling  rates  than  states  localized  at  the  sides  of  the  well.  The 
dependence  of  tunneling  time  on  energy  also  exhibits  steps  related  to  quasi>bound*state  energies  in  a 
complex  way.  Additionally,  modulations  in  the  tunneling  probability  correlate  with  the  motion  of  the 
squeezed  wave  packet  in  the  well.  These  results  have  implications  both  for  the  dynamics  of  squeezed  ra¬ 
diation  fields  and  squeezed  matter  states,  and  examples  of  applications  to  controlled  vibrational-state 
chemistry,  ultralow-noise  measurements  of  weak  signals  with  Josephson  junctions,  and  optical  crystals 
are  discussed. 


PACS  number(s):  03.65.Ge,  32.80.Pj,  42.50.Dv 
L  INTRODUCTION 

■  Theoretical  work  on  the  class  of  minimum-uncertainty 
states  known  as  squeezed  states  was  initiated  as  early  as 
the  1920s  [1],  and  revived  in  the  context  of  quantum  op¬ 
tics  in  the  1960s  [2].  However,  it  as  not  until  1995  that 
experimental  methods  were  developed  to  generate,  propa¬ 
gate,  and  detect  squeezed  states  of  the  radiation  field  [3]. 
This  caused  an  immense  increase  of  interest  in  the  field. 
To  date,  schemes  used  to  generate  squeezed  light  have  in¬ 
cluded  four-wave  mixing  in  atomic  vapor,  four-wave  mix¬ 
ing  in  fibers,  second-order  parametric  interactions  in  non¬ 
linear  crystals,  and  pump-noise-suppressed  lasers  [4-6]. 
It  has  been  shown  that  squeezed  states  are  important  for 
sub-shot-noise-limited  interferometry  and  spectroscopy 
with  sensitivity  below  the  shot-noise  limit,  and  quantum 
nondemolition  measurements.  As  a  consequence  of  these 
applications  and  improved  understanding  of  the  uncer¬ 
tainty  limits  of  quantum  measurement,  interest  has 
grown  recently  in  the  possibility  of  creating  squeezed 
states  of  matter  which  may  similarly  exhibit  unusual  and 
useful  properties  [7]. 

In  this  work,  the  effect  of  squeezing  on  tunneling  from 
a  simple  potential  well  to  a  continuum  of  final  states  is  in¬ 
vestigated.  The  treatment  is  sufficiently  general  to  apply 
to  squeezed  states  of  both  radiation  and  matter  fields. 
However,  it  is  primarily  motivated  by  an  interest  in 
modifications  of  tunneling  phenomena  due  to  localization 
of  matter  wave  packets  within  potential  wells.  Several 
previous  studies  have  analyzed  tunneling  of  squeezed 
states  in  double  well  quantum  oscillators  [8,9],  but  have 
not  been  extended  to  describe  transport.  Here  we  analyze 
the  evolution  of  squeezed  states  in  a  one-dimensional  par¬ 
abolic  well  with  an  infinite  barrier  on  one  side  and  a  finite 
barrier  on  the  other,  explicitly  incorporating  the  concept 
of  transport.  We  show  that  squeezing  of  spatial  coordi¬ 
nates  can  have  a  significant  effect  on  tunneling  probabili- 
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ty  and  introduces  structure  in  the  dependence  of  tunnel¬ 
ing  time  on  energy  of  the  wave  packet  which  is  related  to 
the  occurrence  of  quasibound  states.  Some  of  the  impli¬ 
cations  of  the  sensitivity  of  tunneling  probability  to 
squeezing  in  the  context  of  controlled  vibrational-state 
chemistry,  macroscopic  quantum  tunneling,  and  in  opti¬ 
cal  crystals  are  discussed. 

Calculations  presented  in  the  next  section  use  a  WKB 
procedure  to  obtain  wave  functions  and  the  tunneling 
propagator  for  a  parabolic  well  coupled  to  a  continuum 
region  of  constant  potential  energy.  Analytic  expressions 
are  obtained  for  wave  functions  and  tunneling  probabili¬ 
ties  when  the  system  is  initially  prepared  in  a  squeezed 
state.  Subsequently  we  examine  the  case  of  a  cubic  bar¬ 
rier,  explain  the  relationship  between  localization,  quasi¬ 
bound  states,  and  tunneling  probability,  and  discuss 
several  applications. 

n.  TUNNELING  BY  A  SQUEEZED  STATE 
A.  Wave  fiuctioiis 

The  tunneling  region  considered  in  this  paper  is  depict¬ 
ed  in  Fig.  1,  where  five  distinct  regions  (I-V)  are 
identified  for  convenience.  A  particle  with  average  ener¬ 
gy  E  is  assumed  to  be  initially  localized  in  the  potential 
well  on  the  left  and  tunnels  into  the  classically  allowed  re¬ 
gion  on  the  right  which  is  terminated  by  a  boundary  at 
X  =L.  The  maximum  potential  of  the  barrier  is  denoted 
by  Kq.  The  points  a,  h,  and  c  are  the  classical  turning 
points.  In  the  classically  allowed  region  beyond  the  bar¬ 
rier,  the  potential  falls  to  zero  at  x  =u;,  and  maintains  a 
constant  value  — ^  for  x  >cf,  where  ^  is  assumed  to  be 
positive.  The  potential  energy  in  the  well  region  (region 
II  of  Fig.  1)  is  assumed  to  be  harmonic: 

F|j(x)  =  |m6)oX^  .  (1) 


1050-2947/94/49(l)/32(10)/$06.00 


32 


©1994  The  American  Physical  Society 


49 


TUNNELING  DYNAMICS  OF  SQUEEZED  STATES  IN  A  . . . 


33 


d  L 


FIG.  1.  Schematic  of  the  potential  well  considered  in  this 
work,  coupled  to  continuum  of  states  at  constant  potential  —  ^ 
and  bounded  by  an  infinite  potential  at  x  =L.  Turning  points 
for  energy  E  are  at  x  =a,  b,  and  c. 

Here,  m  is  the  particle  mass,  cdq  is  the  frequency  of  oscil¬ 
lation,  and  X  is  the  particle  position.  A  characteristic 
length  can  be  defined  for  the  linear  harmonic  oscilla¬ 
tor  (LHO)  wave  functions  corresponding  to  this  potential 
weU.  This  localization  length  gives  the  standard  devia¬ 
tion  of  the  particle  position  in  the  LHO  ground  state  and 
is  a  useful  quantity  for  discussing  normalization  and  den¬ 
sity  of  states: 

a  =  \/2ma)o/#  .  (2) 

The  time  evolution  of  a  wave  packet  'F  initially  local¬ 
ized  in  the  well  is  given  by 

'P(x,r)=  r*  (?(x',x;t)'F(x',0)dx'  ,  (3) 

—  00 

where  G(x',x;t)  is  the  propagator  which  accounts  for 
the  evolution  of  the  wave  function  from  its  initial  value 
^(x\0)  at  time  zero  and  position  x'  to  its  later  value 
^(x,r)  at  time  t  and  position  x.  The  propagator 
G{x\x;t)  is  given  by 

k 

Here,  the  sum  is  over  wave  vectors  k  in  the  continuum 
region,  x  >  d,  which  are  defined  by 

k  =  .  (5) 

The  energies  of  the  eigenstates  ^i^lx)  for  the  potential 
F(x)  in  Fig.  1  are  a  discrete  set  derived  from  the  bound¬ 
ary  conditions  —  oo  )=0. 

To  obtain  a  simple  analytic  form  for  the  propagator, 
we  consider  only  states  with  energies  which  lie  below  the 
barrier  peak  Kq  and  use  a  WKB  approach.  We  neglect 
those  states  with  energies  close  to  Vq  for  which  the  WKB 
approximation  with  linear  turning  points  is  not  valid. 
For  example,  the  WKB  approximation  breaks  down  for 
states  within  about  fico^  of  the  barrier  peak  in  a  cubic  po¬ 
tential  containing  roughly  ten  quasibound  states.  For 
moderate  amounts  of  squeezing  and  energies  limited  to 
those  below  Fq,  this  condition  restricts  the  average  ener¬ 
gy  of  the  wave  packet  to  within  two  or  three  below 
the  top  of  the  barrier. 


Upon  applying  the  appropriate  connection  formulas 
for  linear  turning  points  [10],  the  WKB  eigenstates  rf}f^(x) 
with  energies  Q<Ej^<Vq  are  of  the  form 

A(^U^(x)  ,  (6) 

where  i4o  is  a  normalization  constant.  The  unnonnalized 
WKB  solutions  for  the  potential  shown  in  Fig.  1  are 
given  by 

C/^U)=:j^exp[-/J'edx'j,  x<a  (7) 

C/“(jc)=— ^cos  f"gdx'—^  ,  a  <x<b  (8) 

Va  •'a  4 


U^{x)=  ^p=-siti§coi  f^qdx'+-^ 

+  -^cos|cos  f^qdx’  —  ^  ,  c<x<d% 
Vg  [•'c  4  I  f 

In  the  continuum  region  we  find 

if  1  1'^ 

U^(x)=-^  4^008^1+ 

Xcos  k(x—d)+q+<p—^  ,  x>d  (1 

4 

where  the  phase  factor  q>  is  given  by 


^=tan 


r 


-tan|:  . 


The  wave  vectors  for  the  classically  allowed  and  forbid¬ 
den  regions  near  the  ori^  are 

q(x)=aV^[Eic  —  V{x)]/fkt>Q,  x<d  (13) 


Q {x)=a-v/[ K(x)-£fc ]/Aa^  x<d  (14) 

respectively.  In  these  expressions,  we  have  made  con¬ 
venient  use  of  a  barrier  parameter  6  defined  by 


0=exp  [f‘Qdx'  j  . 


The  inverse  of  0^  is  the  tunneUng  probability  for  an  eigen¬ 
state  with  energy  Phases  accumulated  between  turn¬ 
ing  points  in  the  classicaUy  altbwed  region  have  been 
designated  by 

^—f^qdx'  (16) 


T]=f‘‘qdx’  . 
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For  a  parabolic  well,  ^  has  the  particularly  simple  form 


TTE^  _  'TTCOj^ 
^0  ^0 


(18) 


An  examination  of  Eq.  (1 1)  reveals  that  the  probability 
of  the  particle  being  found  in  the  continuum  region  is 
close  to  zero  for  particular  values  of  the  phase  parameter 
I*,  namely  ^={n  +y)7r,  where  n  is  an  integer.  Particles 
with  these  phases  are  strongly  localized  in  the  well  in 
quasibound  states  which,  according  to  Eq.  (18),  occur  al¬ 
most  exactly  at  the  LHO  eigenenergies  and  are  easily 
seen  to  have  narrow  energy  widths  on  the  order  of 
ficoQd~^.  For  energies  below  the  bottom  of  the  well 
<Ei^  <0)y  the  WKB  wave  functions  decay  exponen¬ 
tially  for  X  <c  and  can  be  neglected  in  the  propagator 
- — - - - I 
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B.  Density  of  states  and  the  propagator 

The  determination  of  the  density  of  states  for  evaluat¬ 
ing  the  propagator  requires  similar  considerations.  Ap¬ 
plying  the  boundary  condition  to  Eq.  (11)  im¬ 

plies  A:(L -“d)  +  7;  +  <^-“7r/4  =  (n +|)ir.  The  density  of 
states  is  then  given  by  the  derivative  of  this  quantization 
condition: 


dn 

dk 


TT 


dk  dk 


(21) 


Far  from  resonance,  only  the  first  term  is  significant  in 
the  limit  of  large  L.  However,  as  seen  from  Eq.  (12),  the 
phase  <p  also  displays  strong  resonant  behavior  near  the 
LHO  eigenenergies.  Including  both  of  these  terms  for 
the  general  case  gives  the  density  of  states 


_ 2ka  ^ _ 

dk  IT  40^cos^|4-i0“^sin^|' 


(22) 
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[Eq.  (4)]  for  the  case  of  wave  packets  initially  localized  in 
the  well. 

The  normalization  constant  Aq  is  easily  determined 
from  Eqs.  (6)  and  (11)  for  wave-packet  energies  far  from 
LHO  eigenenergies  when  L  is  large.  However,  near 
quasi-bound-state  resonances,  the  probability  amplitude 
in  the  well  region  is  very  sensitive  to  energy  if  the  condi¬ 
tion 

L»a~^d^  (19) 

cannot  be  satisfied.  For  modest  values  of  10^,  Eq.  (19) 
implies  that  L  must  be  on  the  order  of  centimeters  if  the 
localization  of  the  wave  packet  within  a  well  a  few 
angstroms  in  size  is  to  be  meaningful.  As  a  consequence, 
we  consider  oflf-resonant  and  near-resonant  cases  sepa¬ 
rately  (Appendix  A),  finding 


«+}) 

I'asTrin  -l- j)  . 


(20) 


The  propagator  can  now  be  evaluated  by  replacing  the 
summation  over  k  in  Eq.  (4)  by  an  integral  over  the  densi¬ 
ty  of  states.  Applying  the  results  of  Eqs.  (6),  (20),  and 
(22)  yields  the  following  expression  for  the  propagator, 
valid  both  for  on-resonance  and  off-resonance  conditions: 


ke  '^*'^V**(x')£/*(x) 
40^cos^^ + 


x>d  . 


(23) 

The  denominator  in  the  integrand  of  Eq.  (23)  becomes 
strongly  resonant  at  quasi-bound-state  energies.  Only 
states  near  these  energies  make  important  contri¬ 

butions  to  the  integral.  Moreover,  since  the  initial  wave 
functions  are  assumed  to  be  strongly  localized  in  the  well, 
these  states  can  be  replaced  by  the  LHO  Hermite- 
Gaussian  wave  functions  with  appropriate  normalization 
(Appendix  A).  With  a  change  of  variable  from  k  to  co, 
the  propagator  becomes 


G(x\x;t)  = 


a 


2  {-i)"^fr;r”°(x') 


( n  + 1  )<an 


-1/2^ 


'0  n  =0 


^i[k{x  —d)-^7f—iir/4] 

- - - hc.c. 

[2dcos^—{i/2)d  'sinl" 


da  , 

x'<c,  x>d 


(24) 


Here,  the  remaining  summation  is  over  the  N  quasibound 
states  that  lie  below  the  barrier  for  which  the  WKB  and 
parabolic  approximations  are  valid. 

The  denominators  in  the  integrand  of  Eq.  (24)  can  be 
expanded  about  each  resonance  and  are  well  approximat- 
^  by  simple  poles.  For  ^«(«  -H  j)ir,  the  denominators 
can  be  written  as 


1 


20cos|-±(//2)0-'sin| 


s( -!)"  +  > _  "  " _ 

[ca— (« -l-i)6)o]±/(r,/2) 


where 
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r„  = 


6>o 

■ 


(26) 


In  this  approximation,  we  ignore  frequency  shifts  on  the 
order  of  o)Qd~*  and  higher-order  imaginary  contributions 
which  introduce  corrections  to  these  first-order  quasi¬ 
bound-state  decay  rates.  Note  that  we  have  also  made 
explicit  use  of  the  parabolic  well  assumption. 

To  complete  the  evaluation  of  Eq.  (24),  we  expand  the 
phase  7]  about  frequency  Q}=(n  +i)<uo.  From  Eq.  (17), 
we  find 


+7)0)0]  , 


(27) 


classically  allowed  region  between  the  barrier  and  the 
continuum;  t„  is  the  respective  classical  transit  time 
across  the  same  region. 

Similarly,  the  wave  vector  k  may  be  expanded  about  a, 
yielding 


*=*(.+r«[o)-(«+{)<i)o] . 
where 

]l/2 


k„=a 


,  1  .  (f> 
2 


and 


(30) 


(31) 


where 

Tt„=a  /V(«+|)-nx')/Au(^^ '  (28) 

and 


T.= 


a 

2mn 


f. 


dx' 


c  'v/(n+i)-F(x')/Auo 


(29) 


Here,  ij„  is  the  phase  accumulated  from  traversing  the 


(32) 


By  extending  the  limits  of  integration  in  Eq.  (24)  to 
infinity  and  making  use  of  these  results,  the  integral  can 
easily  be  evaluated  by  the  method  of  residues.  Retaining 
only  the  term  representing  a  wave  traveling  to  the  right, 
and  removing  the  slowly  varying  term  k  from  the  in¬ 
tegral,  the  final  expression  for  the  propagator  is 


G(x',x;t)= 


a 

-1/2  ^ 

n  -(-  —  -I — ^ 

Air 

n=0 

2  fioio 

-1/4 


Q[t-r„-r„{x  -d)\ 


(33) 


Here,  0(x)  is  the  Heaviside  step  function.  The  arbitrary 
exponential  phase  factor  tt/A  in  Eq.  (24)  has  been  ig¬ 
nored. 

The  propagator  in  Eq.  (33)  consists  of  a  superposition 
of  exponentially  decaying  traveling  waves  which  are  de¬ 
layed  by  corresponding  classical  transit  times  from  the 
edge  of  the  barrier  to  the  continuum  region.  Each  of  the 
quasibound  states  decays  with  the  usual  WKB  decay  rate 
given  by  Eq.  (26).  The  sharp  tum-on  of  each  quasibound 
state  is  related  to  the  simple  pole  approximation  used  to 
evaluate  the  propagator.  We  point  out  that  these  results 
are  also  valid  when  the  continuum  lies  above  the  bottom 
of  the  weU  (^  <0)  if  the  lower  limit  in  the  summation  in 
Eq.  (33)  is  replaced  by  the  lowest  LHO  level  which  lies 
above  the  continuum  energy.  True  LHO  bound  states  in 
the  well  which  lie  below  the  continuum  decay  exponen- 
tiaUy  in  the  region  x  >  6  in  Fig.  1  and  make  negligible 
contributions  to  the  probability  beyond  the  barrier. 


C.  Tnnneling  probability 

Application  of  Eqs.  (3)  and  (33)  to  an  arbitrary  wave 
packet  composed  of  LHO  eigenstates  with  probability 
amplitudes  reveals  that,  beyond  the  barrier,  'Kx,/)  in 
Eq.  (33)  is  simply  a  superposition  of  tunneling  com¬ 
ponents  weighted  by  their  appropriate  c„.  The  probabili¬ 
ty  of  escape  for  a  wave  packet  initially  localized  in  the 
well  is  defined  to  be 


PU)=  r\'*(x',t)\^dx'  .  (34) 

d 

Ignoring  the  rapid  interference  modulations  between 
different  quasibound  states,  this  gives  a  tunneling  proba¬ 
bility  of 

P(t)  =  — - - - ,  (35) 

when  normalized  to  the  probability  of  the  initial  wave 
packet  lying  below  the  barrier.  This  result  is  easily  inter¬ 
preted.  The  tunneling  probability  is  simply  the  sum  of 
the  individual  probabilities  for  eacjti  of  the  exponentially 
decaying  quasibound  states  to  tunnel  through  the  barrier, 
weighted  by  the  probabUity  of  the  wave  packet  to  be  in 
that  particular  quasibound  state.  As  t  — oo ,  the  probabil¬ 
ity  in  Eq.  (3S)  approaches  unity. 

m.  TUNNELING  OF  SQUEEZED  STATES 
A  Tniuielhig  tines 

To  investigate  the  significance  of  squeezing  on  tunnel¬ 
ing  times,  we  calculate  the  tunneling  probability  in  Eq. 
(35)  for  the  case  in  which  the  initial  state  4'(x,0)  is  a 
squeezed  state.  To  do  this,  we  evaluate  d  and  tj  for  the 
specific  case  of  a  cubic  barrier  coupled  to  a  flat  continu- 
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um  (Fig.  1).  The  barrier  is  assumed  to  have  a  height 
Fo  =  10Ayo»  and  we  take  the  number  of  corresponding 
quasibound  states  to  be  iV  =  10.  Details  of  the  evaluation 
of  barrier  quantities  are  presented  in  Appendix  B. 

We  consider  only  pure  amplitude  and  phase-squeezed 
states.  The  initial  squeezed  wave  packet  can  be 
represented  as  [1 1] 

.  (36) 

The  standard  deviation  Ax  of  wave-packet  position 
within  the  parabolic  well  is  (a5)”*,  where  S  is  related  to 
a  squeezing  parameter  R  which  is  a  linear  measure  of 
squeezing  and  defined  by  the  relation  [11].  If 

5  <  1,  the  wave  packet  is  phase  squeezed  and  character¬ 
ized  by  a  Ax  greater  than  the  ground-state  standard  devi¬ 
ation.  If  S>  1,  amplitude  squeezing  applies,  and  Ax  is 
smaller  than  the  ground-state  standard  deviation.  The 
initial  mean  displacement  in  2p{aS)~^, 

The  initial  squeezed  state  may  easily  be  represented  in 
terms  of  LHO  states  with  probability  amplitudes  which 
can  be  calculated  by  various  methods  [11,12].  Tunneling 
escape  times  r^,  arbitrarily  defined  by  the  condition 

P(r,  )=!-«-',  (37) 

may  then  be  calculated  by  direct  application  of  Eq.  (35). 
This  definition  of  escape  time  permits  direct  comparison 
with  the  purely  exponential  decay  of  LHO  eigenstates. 
The  tunneling  escape  times  for  various  degrees  of  squeez¬ 
ing  versus  average  energy  are  plotted  in  Fig.  2,  together 
with  tunneling  times  for  eigenstates  of  the  parabolic  well. 
Figure  3  reveals  that  the  tunneling  times  for  squeezed 
states  reach  a  minimum  at  5  « 1.5  for  an  average  energy 
£=4.5^0.  Interestingly,  this  does  not  correspond  to 
the  case  of  a  pure  coherent  state  (S  =  1.0).  Instead,  it 
corresponds  closely  to  the  value  of  S  at  which  energy 
fluctuations  reach  a  minimum,  as  indicated  in  Fig.  4  and 
discussed  below. 

To  understand  first  how  squeezing  affects  the  energy 
distribution,  we  consider  the  relationship  between  squeez¬ 
ing  and  the  energy  variance  <(A^)^).  This  is  given  by 


FIG.  2.  Tunneling  escape  times  versus  energy  for  cubic  bar¬ 
rier,  plotted  for  various  values  of  squeezing:  5  =0.5  (solid),  1.0 
(dashed),  1.5  (dotted),  and  3.0  (chain).  Barrier  and  continuum 
parameters  are  Vq  =  lOfifu©,  N  =  10,  and  0  =  2^o- 


FIG.  3.  Example  of  tunneling  escape  times  versus  squeezing 
parameter  for  an  arbitrary  fixed  energy  E  =4.  SfuoQ. 

the  expression  [12] 

<(A£)2>  =  j(«£ao)2[(8)S^+l)S-‘+S‘‘-2]  .  (38) 

The  energy  of  a  squeezed  state  is  given  by  a  similar  ex¬ 
pression, 

<£>  =  iftao[(4^^+l)S“^+S^]  .  (39) 

The  average  energy  in  Eq.  (39)  depends  both  on  the  initial 
mean  displacement  of  the  squeezed  state  2p{aS)~^  as 
well  as  the  position  and  momentum  uncertainties,  Ax  and 
Ap.  These  in  turn  are  proportional  to  5“^  and  5,  respec¬ 
tively.  These  uncertainties  contribute  additional  kinetic 
and  potential  energy  to  the  wave  packet.  Increasing  its 
energy  above  that  of  a  pure  coherent  state  with  displace¬ 
ment  0, 

From  Eq.  (39),  we  see  that  for  small  displacements  0 
the  energy  is  dominated  by  uncertainties  in  position  and 
momentum.  For  small  energies  (2iSS'"*«0),  these  uncer¬ 
tainties  contribute  equally  to  potential  energy  ( Ax  )^  and 
kinetic  energy  (Ap)^.  Thus,  for  small  energies,  the  total 
fluctuation  in  energy  is  minimized  when  potential  and 
kinetic  energy  fluctuations  are  simultaneously  minimized, 


FIG.  4.  Energy  *  fluctuations  of  a  squeezed  state  versus 
squeezing  parameter  S.  The  three  curves  correspond  to 
different  particle  energies:  E  =2.0  (solid),  5.0  (dashed».  and  50 
(dotted). 
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or  S  is  close  to  one.  On  the  other  hand,  for  larger  initial 
displacements  (2^5~*>0),  the  potential  energy  acquires 
contributions  not  only  from  position  and  momentum  un¬ 
certainties,  but  from  the  mean  displacement  itself.  Fur¬ 
thermore,  large  wave-packet  displacements  increase  the 
slope  of  the  potential  energy  curve,  amplifying  the  contri¬ 
butions  of  positional  fluctuations  Aa:.  TTius,  at  large  aver¬ 
age  energies  {27jS~^»0\  fluctuations  of  potential  ener¬ 
gy  are  weighted  more  heavily  than  fluctuations  of  kinetic 
energy.  In  this  way,  energy  fluctuations  can  be  expected 
to  minimize  for  a  value  of  S  greater  than  one,  in  agree¬ 
ment  with  an  algebraic  minimization  of  Eq.  (38)  with 
respect  to  S  which  furnishes  the  result  1  + 

The  corresponding  minimization  of  tunneling  times 
occurs  as  a  result  of  the  LHO  occupation  probability  dis¬ 
tributions  narrowing  to  a  small  number  of  LHO  quasi¬ 
bound  states  as  S  approaches  5^  (Fig.  5).  Quasi- 
bound-state  tunneling  probabilities  for  neighboring  levels 
vary  by  several  orders  of  magnitude,  so  a  wave  packet 
composed  of  a  wide  distribution  of  LHO  states  requires 
many  very  long-lived  quasibound  states  to  decay  before 
the  tunneling  probability  becomes  appreciable.  However, 
for  squeezed  wave  packets  with  only  a  few 

short-lived,  quasibound  states  need  to  decay  before  the 
wave  packet  has  tunneled  significantly.  This  also  ex¬ 
plains  why  LHO  eigenstates  have  the  shortest  tunneling 
escape  times  for  a  given  energy,  as  seen  in  Fig.  2. 

In  addition  to  the  sensitivity  of  wave-packet  escape 
times  on  squeezing,  the  tunneling  escape  times  plotted  in 
Fig.  2  also  exhibit  interesting  stepUke  features  near 
quasi-bound-state  energies.  These  steps  are  related  to  the 
decomposition  of  the  initial  squeezed  state  into  a  super¬ 
position  of  discrete  quasibound  states.  A  squeezed  state 
dominated  by  a  single  quasibound  state  exhibits  tunneling 
behavior  dominated  by  that  level.  If  the  average  energy 
is  increased  by  some  fraction  of  then  additional  lev¬ 
els  must  be  occupied  significantly  to  account  for  this  ad¬ 
ditional  energy.  Because  of  the  very  large  differences  in 
tunneling  rates  between  neighboring  LHO  levels, 
significant  changes  in  tunneling  rates  are  to  be  expected 
as  higher  LHO  levels  become  significantly  occupied. 


FIG.  5.  Probability  distribution  for  squeezed  states 
versus  occupation  number  n,  plotted  for  various  values  of 
squeezing:  5=0.5  (circles),  1.5  (squares),  and  3.0  (triangles) 
with  arbitrary  fixed  energy  E  =4.5fut)Q,  The  curves  are  to  guide 
the  eye. 


For  large  amounts  of  phase  squeezing  (S  =0.5,  Fig.  2) 
the  steps  in  tunneling  time  become  less  prominent  and  do 
not  occur  at  every  LHO  eigenenergy.  Under  these  condi¬ 
tions,  the  broad  energy  distributions  of  phase-squeezed 
states  cover  more  than  one  LHO  level  spacing  so  that  the 
correspondence  between  particular  quasibound  states  and 
individual  steps  is  lost.  The  tunneling  time  also  exhibits 
steplike  features  as  a  function  of  squeezing  for  fixed  ener¬ 
gies,  as  depicted  in  Fig.  3.  These  features  similarly  arise 
from  changes  in  occupation  of  LHO  discrete  levels  as  the 
amount  of  squeezing  is  varied. 

B.  Coherences 

Coherent  effects  become  evident  when  the  properties  of 
the  timneled  wave  function  can  be  resolved  on  time  scales 
on  the  order  of  and  length  scales  on  the  order  of 
For  example,  the  probability  \tfix,t)\^  for  a  pure 
coherent  state  develops  deep  modulations  at  the  wave- 
packet  oscillation  frequency  oq  as  shown  in  Fig.  6.  These 
modulations  arise  from  interferences  between  adjacent 
quasibound  states. 

Wave-packet  squeezing  can  emphasize  higher  harmon¬ 
ics  of  6>o  and  alter  the  shape  of  tunneling  modulations. 
These  effects  are  apparent  in  Rgs.  7{a)-7(d),  aiiere 
is  plotted  for  x=d  with  various  degrees  of 
squeezing  at  a  fixed  energy  close  to  the  LHO  ground-state 
energy.  Figure  7(a)  depicts  a  pure  coherent  state  with 
modulations  at  Wq,  consistent  with  the  purely  oscillatory 
motion  of  the  wave-packet  mean  <x(f)>  in  the  well.  As 
the  wave  packet  becomes  amplitude  squeezed,  the  wave- 
packet  variance  < [Ax(t)]*)  begins  to  oscillate  at  Ioq,  so 
that  second  harmonic  modulations  begin  to  appear  in  the 
tunneling  probability  beyond  the  barrier  [Fig.  7(b)].  Simi¬ 
larly,  phase  squeezing  introduces  second  harmonics,  but 
with  the  phase  reversal  evident  in  Hg.  7(c). 

With  maximum  squeezing  (P=0),  wave  packets  local¬ 
ized  in  the  well  consist  of  superpositions  of  even  harmon¬ 
ics  only,  with  the  result  that  breathing-type  oscillations  in 
([Ax{t)f)  occur  purely  at  2wq  as  shown  in  Fig.  7(d).  At 
higher  energies,  these  coherent  effects  become  less  pro- 
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FIG.  6.  Probability  |^d,t)p  of  particle  detection  at  point  d 
beyond  the  barrier  as  a  function  of  time  for  a  coherent  state 
with  5  =  1.0  and  energy  E=0.6fiG>o.  Barrier  and  continuum 
parameters  are  and  ^=2fifiio* 
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and  initial  displacements  p  for  particle  energies  exceeding  the 
ground-state  energy  Eg=0.5fio}Q  by  only  =  lO~^fut)Q:  (a) 
S  =  1,0000,  i9=10-^;  (b)  5  =  1.0001,  ^=9.95 XIO'^  (c) 
5  =  0,9999,  )9=9.95X10‘^  and  (d)  5  =  1.0010,  i5=0.  Barrier 
and  continuum  parameters  are  Vq  =4^o»  iV  =4,  and  <l>=2fuoQ. 


nounced  as  excited  quasibound  states  become  progres¬ 
sively  more  populated.  This  occurs  because  d~^  increases 
exponentially  with  energy  level  n  and  dominates  the 
behavior  of  the  tunneling  probability  contributions 
at  higher  energies,  reducing  modulation  con¬ 
trast.  Modulations  become  large  only  when  the  probabil¬ 
ities  decrease  rapidly  enough  to  counterbalance  the 
d~^  term  in  For  squeezed  states  tunneling 

through  a  cubic  barrier,  this  occurs  at  average  particle 
energies  approaching  that  of  the  ground  state, 

IV.  DISCUSSION  AND  CONCLUSIONS 

The  results  obtained  in  this  paper  have  interesting  im¬ 
plications  for  quantum  control  [13]  of  chemical  systems 
with  eflfective  potentials  of  the  general  form  shown  in  Fig. 
1,  where  an  energy  barrier  separates  the  initial  from  the 
final  state.  For  example,  simple  molecules  such  as  dia¬ 
tomic  alkali  metals  occasionally  have  excited-state  poten¬ 
tials  with  centrifugal  barriers  separating  shallow  quasi¬ 
bound  regions  from  dissociative  continua.  The  state 
of  Cs2  accessed  at  766.7  nm,  for  example,  has  a  potential 
of  this  type  [14].  In  principle,  for  excited-state  potential 
wells  deep  enough  to  contain  several  vibrational  states, 
dissociative  dynamics  would  then  depend  on  the  detailed 
nature  of  the  initial  wave  packet. 

Our  results  are  most  relevant,  however,  to  cases  of  very 
deep  potential  wells  with  part  or  all  of  the  well  located 
above  the  dissociative  asymptote.  This  situation  is  more 
typical  of  predissociative  states  of  simple  molecules  like 
IBr  (electronic  predissociation  [15])  and  HgH  (vibrational 
dissociation  [16]).  Examples  among  the  diatomic  alkali 
metals  include  the  C  state  of  Cs2  [17]  and  the  6  ^2^ 
state  of  Na2  [18]  which  form  adiabats  through  avoided 
crossings  with  repulsive  states  similar  to  the  model  con¬ 
sidered  here.  The  influence  of  tunneling  was  specifically 
considered  as  early  as  the  work  of  Farkas  and  Levy  [19] 
on  AlH.  In  such  systems,  our  calculations  indicate  that 
initial-state  preparation  strongly  affects  the  evolution  of 
excited-state  reactions.  Both  the  amplitude  and  phase  of 
individual  quasibound  states  encompassed  in  the  initial 
superposition  state  affect  the  rate  of  tunneling  through 
the  barrier.  Reaction  time  scale  and  yield  should  there¬ 
fore  be  amenable  to  manipulation  using  incident  light 
pulses  of  durations  less  than  a  vibrational  period  with 
tailored  amplitude,  bandwidth,  central  frequency,  and 
chirp  to  slow  down  or  speed  up  chemical  dynamics.  In 
optically  initiated  unimolecular  reactions,  ultrafast  dissi¬ 
pation  of  excitation  to  intramolecular  vibrational  reser¬ 
voirs  should  be  impeded  by  appropriate  phase  and  chirp 
manipulation  of  incident  puls^.  This  suggests  useful  ex¬ 
tensions  of  the  experimental  techniques  developed  recent¬ 
ly  [20]  for  molecular  wave-packet  preparation,  to  exploit 
optical  phasing  in  the  creation  of  squeezed  matter  states 
analogous  to  squeezed  states  of  the  radiation  field  [11] 
with  correspondingly  modified,  nonclassical  dynamics. 

The  modification  of  tunneling  dynamics  by  squeezing 
also  has  interesting  consequences  for  dc-biased  Josephson 
junctions  and  Josephson  interferometers  and  is  particu¬ 
larly  germane  to  precision  measurements  of  magnetic  flux 
or  voltage  close  to  the  shot-noise  limit  with  these  struc- 
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tures  [21].  Here,  the  relevant  variable  is  the  Cooper  pair 
center-of-mass  phase  difference  8  in  the  BCS  ground-state 
wave  functions  on  either  side  of  the  Josephson  barrier 
[22],  The  dynamics  of  the  phase  difference  8  arise  from  a 
quantum-mechanical  “binding  energy”  F(8)  due  to  cou¬ 
pling  of  Cooper  pairs  on  either  side  of  the  barrier  by  the 
tunneling  interaction  [23].  These  dynamics  have  been 
studied  extensively  in  both  the  classical  regime  [24]  and 
more  recently  in  the  quantum  regime  at  low  temperatures 
where  quantized  energy  levels  observed  in  interactions 
with  microwave  radiation  [25]  and  tunneling  of  the  phase 
difference  8  from  potential  minima  in  K(8)  have  been  ob¬ 
served  [26]. 

The  tunneling  dynamics  of  the  phase  difference  8  and 
corresponding  junction  voltage  V=(ii/le)b  across  a 
Josephson  junction  with  capacitance  C  can  be  described 
by  an  effective  quantum-mechanical  Hamiltonian  consist¬ 
ing  of  a  potential  energy 

K(8)=  ~Ey(cos8+8///y)  (40) 

and  a  kinetic  charging  energy  (len  )^/2C,  where  2en  is  the 
amount  of  charge  transferred  across  the  junction, 
Ij~2eEj /fly  and  £}  is  a  constant  related  to  the  tunneling 
interaction  matrix  elements  and  density  of  states 
[23,27,28].  When  the  applied  current  I  is  less  than  Ij  in 
Eq.  (40),  the  potential  F(8)  is  locally  well  approximated 
by  a  cubic  potential  of  the  form  shown  in  Fig.  1  with 
<^»ikoQ.  The  preparation  of  squeezed  states  tffib)  should 
then  be  feasible  by  application  of  a  constant  bias  current 
and  a  microwave  current  at  2o}q  to  modulate  the 
resonant  frequency  of  the  “quantum  well”  [25] 
Q)Q==^(2V^eIjC/fi)^^{\--I/Ij)^^^  together  with  seed 
currents  at  o)q  to  provide  additional  control  of  mean 
wave-packet  displacement  and  phasing.  This  is  closely 
analogous  to  optical  [5]  and  microwave  squeezing  [29]  by 
parametric  interactions. 

The  results  of  this  work  indicate  that  large  amounts  of 
phase  squeezing  (5  » 1 )  would  strongly  inhibit  tunneling 
decay  in  a  potential  of  the  form  given  in  Eq.  (40),  and 
permit  improved  ultralow-noise  measurements  of  weak 
signals  at  low  temperatures  by  squeezed  states  whose 
noise  properties  are  well  known  to  be  highly  affected  by 
loss  [30].  Furthermore,  we  expect  these  results  to  be  im¬ 
portant  in  fundamental  dynamics  of  coupled  junctions 
where  wave-packet  loca^tion  can  be  expected  to 
enhance  or  inhibit  tunneling  in  a  fashion  reminiscent  of 
numerical  simulations  of  forced  double  potential  wells 
[31]. 

Optical  crystals,  formed  by  trapping  laser-cooled 
atoms  in  periodic  optical  potential  wells,  constitute 
another  system  in  which  squeezing  of  atomic  states  could 
significantly  affect  tuimeling.  In  this  case,  tunneling 
occurs  as  nonclassical  motion  of  escaping  atoms  or  those 
moving  between  minima  of  the  optical  field  configuration. 
Recent  results  have  demonstrated  quantized  motion  of 
atoms  trapped  within  optical  potential  wells,  indicated  by 
the  presence  of  vibrational  sidebands  in  resonance 
fluorescence  spectra  or  in  stimulated  Raman  spectra  [32]. 

J ust  as  for  quasibound  states  of  molecules,  the  induced  vi¬ 
brational  manifolds  of  such  atoms  should  be  amenable  to 


preparation  of  squeezed  vibrational  states  [7].  The  re¬ 
sults  of  the  present  work  indicate  that  squeezing  of 
trapped  atoms  could  strongly  influence  their  tunneling 
rates  into  neighboring  wells  as  well  as  their  escape  from 
the  trap. 

In  summary,  we  have  calculated  tunneling  times  for 
squeezing  states  initially  confined  to  a  harmonic  potential 
well  coupled  to  a  continuum.  The  fastest  tunneling  rates 
are  observed  for  wave  packets  prepared  with  smaU 
amounts  of  amplitude  squeezing,  but  their  rates  never 
exceed  those  of  pure  LHO  eigenstates.  On  the  other 
hand,  tunneling  is  inhibited  by  large  squeezing.  We  have 
also  elucidated  the  relationship  between  the  energy 
dispersion  of  wave  packets  and  tunneling  rates,  as  well  as 
more  subtle  steplike  behavior  arising  from  the  discrete 
nature  of  quasibound  states  in  the  well.  Oscillatory 
behavior  of  the  wave  packet  can  give  rise  to  modulations 
in  the  tunneling  wave-function  probability  which  are 
most  visible  when  relative  contributions  to  excited-state 
populations  decrease  exponentially  with  increasing  vibra¬ 
tional  quantum  number.  Finally,  applications  of 
squeezed  matter  wave  packets  have  been  identified  for 
control  of  excited-state  chemical  reactions  by  appri^ri- 
ately  tailored  ultrafast  light  pulses,  modification  of  mtu:- 
roscopic  tunneling  dynamics  in  Josephson  junctions,  ind 
inhibition  of  tunneling  of  trapped  atoms  in  optical  crys¬ 
tals. 
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APPENDIX  A 


The  question  of  normalization  of  wave  functions  may 
be  handled  by  first  matching  the  WKB  solutions  (x) 

to  the  normalized  LHO  eigenstates  in  the  semiclassical 
limit.  By  explicit  integration,  the  WKB  solution  in  Eq. 
(8)  for  a  parabolic  well  yields 


(z)=  (2a)'^ 
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where  z=ax.  The  phase  d  is  given  by 
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In  the  limit  of  large  quantum  numbers,  this  result 
matches  the  asymptotic  expansion  for  the  parabolic 
cylinder  functions  DJ—z)  to  within  a  constant  of  pro- 
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portionality  [33].  For  integer  values  of  n,  the  parabolic 
cylinder  functions  reduce  to  Hermite-Gaussian  LHO 
solutions.  Comparing  the  WKB  solution  in  Eqs.  (Al) 
and  (A2)  to  the  normalized  asymptotic  LHO  solutions  for 
integer  n  reveals 

f/,U)  =  2(-l)V/^a-'tfrJ:«0(;c),  ;t<c  (A3) 

where 

^i.HO( ^ ^ 2 - ^ /2( 2ir ) “ ' ^*a ' ~ /Vl )  . 

(A4) 

Here,  H„(x)  are  Hermite  polynomials  of  order  n. 

We  can  now  calculate  the  normalization  constant  Aq 
for  the  on-resonance  and  oflf-resonance  cases  by  integrat¬ 
ing  in  Eq.  (6)  over  all  space.  When  the  resonance 

condition  |'=ir(n  +  y)  is  not  met,  only  the  wave-function 
amplitude  in  continuum  region  contributes  significantly 
to  the  total  wave-function  probabihty.  Integrating 
\Uk(x)\^  over  the  continuum  region  gives  the  off- 
resonant  result  in  Eq.  (20).  However,  close  to  quasi¬ 
bound-state  resonances,  the  wave-function  ampUtude  in 
the  well  region  becomes  appreciable.  Replacing  the 
WKB  wave  functions  in  the  well  region  [Eqs.  (7)-(9)] 
with  the  LHO  wave  functions  normalized  to  the  WKB 
solutions  in  Eq.  (A3),  we  find  i*  approximately 

Aq=  4ira~^f^ 

—  oo 

-1/2 

+  J^\Ul(x')\^dx'  ,  |«ir(/i+i).  (A5) 

u  “  ^ 

The  small  contributions  just  beyond  the  barrier  (region 
IV)  have  been  neglected.  Evaluating  Eq.  (A5)  explicitly 
yields  the  result  for  the  on-resonant  case  in  Eq.  (20). 

APPENDIX  B 

In  this  appendix,  we  compute  the  barrier  parameters 
for  the  specific  case  of  a  cubic  potential.  The  cubic  po¬ 
tential  in  Fig.  1  has  the  general  form 

V { x) = (27 Vq / A)xHw —x)/w^  .  (Bl) 


At  the  LHO  eigenenergies  E  =fuoo{n  -bi),  the  tunneling 
parameter  [Eq.  (15)]  is 

0„=exp  af‘-\/V(x)/imo-(n  +±)dx’  ,  (B2) 

where  b  and  c  are  the  barrier  turning  points.  To  obtain 
accurate  results  for  intermediate  energies,  the  integral  in 
(B2)  is  evaluated  exactly.  This  can  be  performed  by  fac¬ 
toring  the  cubic  polynomial  in  the  radical  of  the  in¬ 
tegrand  in  Eq.  (B2).  The  zeros  of  the  polynomial  are  sim¬ 
ply  the  three  turning  points  a,  b,  and  c.  The  integral  can 
be  expressed  in  terms  of  the  Gauss-hypergeometric  func¬ 
tion  [34].  The  resulting  tunneling  parameter  becomes 


d„=exp(/,),  (B3) 

where  I„  is 


y  „  (  13,.  P3-3’2 

[  2  2  y^-yx  ) 

(B4) 

The  parameters  J'3>>’2>yi  are  the  turning  points  in 
units  of  w  and  are  roots  of  the  cubic  equation 

B^yHy  —  l}+{n  -l-i)=0  . 

(B5) 

The  constant  B  is  given  by 

B  =3/27Ko/4ft»o  • 

(B6) 

The  phase  parameter  17,  and  the  classical  delay  time  t„, 
given  by  Eqs.  (28)  and  (29),  are  most  easily  evaluated  by 
assuming  the  potential  in  region  IV  of  Fig.  1  is  linear. 
We  simply  state  the  results: 

[  3B  J  p  2  Auo 

(B7) 

and 

\  1  A  V'* 

r„=2Bmo  n+}  +  -^  .  . 

2  ft»o 
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ABSTRACT 

We  present  results  on  the  first  application  of  upconversion  dynamics  for  short  pulse  generation  in 
visible  solid  state  lasers.  We  also  describe  prospects  for  passive  mode-locking  mechanisms  based  on 
cooperative  nonlinearities  and  instabilities. 


I.  INTRODUCTION 

In  recent  years,  Coulombic  interactions  between  impurity  atoms  in  concentrated  laser  media  have  - 
attracted  attention  in  the  context  of  continuous-wave  upconversion  lasers  and  enhanced  operation  of 
conventional  solid  state  lasers  [1,2].  However,  at  least  two  aspects  of  cooperative  dynamics  of  rare 
earth  ions  in  concentrated  solids  have  interesting  implications  for  ultrashort  pulse  generation.  In  this 
paper,  the  unique  ability  of  cooperative  upconversion  to  furnish  visible  sources  of  mode-locked  pulses 
in  solids  is  first  outlined.  Then,  the  potential  use  of  cooperative  nonlinearities  as  a  new  mechanism  for 
passive  mode-locking  of  Nd  and  other  solid  state  lasers  is  explored. 

The  development  of  continuous-wave  and  cw  mode-locked  sources  at  short  wavelengths  in  solid  state 
media  is  a  challenging  problem.  Many  factors  contribute  to  this,  but  chief  among  them  are  losses 
associated  with  transient  color  center  and  band  edge  absorptions  at  pump  and  lasing  wavelengths. 

This  problem  is  compounded  by  the  scarcity  of  pump  sources  at  short  wavelengths.  One  approach 
which  overcomes  such  difficulties  is  the  use  of  upconversion,  shown  in  recent  years  to  permit  short 
wave  laser  operation  sustained  by  long  wave  excitation.  Cooperative  upconversion  in  particular  is 
considered  here  as  a  means  of  furnishing  useful  gain  in  the  visible  and  ultraviolet  regions  for  direct 
generation  of  short  pulses  by  active  mode-locking  [3]. 

Additionally,  we  explore  the  potential  of  the  nonlinear  susceptibility  generated  by  cooperative 
processes  as  a  novel  mechanism  for  passive  mode-locking.  Both  die  nonlinear  response  due  to 
population  dynamics  and  that  due  to  delocalization  of  rare  earth  intwactions  are  examined  for  prospects 
of  new  sources  with  self-starting,  passively  mode-locked  operation  achieved  through  cooperative 
dynamics. 

EXPERIMENTS 

Using  apparatus  ^scribed  frcviously  [4],  cw  mode-locked  pulse  generation  in  the  green  spectral 
region  was  demonstrated  using  cooperative  upconversion  of  Erbium  ions  in  5%  EnLiYF4  excited  by 
infrared  radiation  at  1.5  nm.  Infiar^  pumping  was  provided  by  a  cw  NaQ  color  center  laser  tuned 
with  a  grating  to  the  first  absorption  resonance  of  trivalent  Erbium.  A  3  mm  thick  laser  crystal  was 
cooled  to  9K  inside  a  dewar  which  also  contained  all  of  the  cavity  optics  except  the  output  coupler  in  a 
standng  wave  geometry  with  astigmatism  correction.  A  self-resonant,  piezo-acoustic  mode-locker 
fabricated  from  LiNb03  was  inserted  at  Brewster's  angle  in  the  collimated  arm  of  the  cavity  external  to 
the  dewar.  Using  active  AM  mode-locking  at  a  repetition  rate  of  1 19  MHz,  pulse  trains  were  then 
generated  continuously  at  544  nm  when  the  cavity  length  was  adjusted  to  within  50  |im  of  the 


synchronous  length  near  60  cm.  Instrument- limited  pulsewidths  of  180  ps  were  observed  in  real  time 
using  a  fast  photodiode  and  a  20  GHz  sampling  oscilloscope. 

In  a  second  set  of  e.xperiments,  light  from  a  cw  Nd;  YAG  laser  was  focused  into  a  single  mode  silica 
fiber  doped  with  27c  Tm3+.  A  maximum  of  70  mW  at  1 .06  (im  was  coupled  into  the  fiber.  The 
fundamental  Nd  wavelength  is  not  in  resonance  with  any  ground  state  absorption  of  Tm^"^  but 
nevertheless  caused  strong  excitation,  visible  to  the  naked  eye.  Blue  upconversion  emission  from 
excited  Thulium  ions  was  monitored  through  interference  filters  at  the  exit  of  a  10  m  fiber,  and 
measurements  made  to  determine  the  dependence  of  upconversion  on  input  power.  At  the  same  time, 
evidence  was  sought  for  a  measurable  threshold  in  the  upconversion  emission,  indicative  of  avalanche 
upconversion  [5,6]  rather  than  step-wise,  multi-photon  absorption  among  excited  states  [7]. 

RESULTS 

Cooperative  upconversion  among  Er  ions  in  LiYF4  was  able  to  sustain  continuous-wave,  mode-locked 
operation  of  a  solid  state  laser  in  the  green  spectral  region  at  liquid  helium  temperatures.  The  output 
pulse  train  is  shown  in  Fig.  1. 
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Figure  1.  Mode-locked  output  versus  time  at  0.5440  pm  for  an  infrared  pump  power  at  1.5  pm  of  500 
mW  power  focused  to  a  spot  radius  of  18  pm  in  the  gain  medium.  T=9K. 

The  dependence  of  pulsewidth  on  output  mirror  position  is  shown  in  Fig.  2.  An  interesting  result  of 
restricting  the  inversion  mechanism  to  cooperative  processes  was  that  the  tendency  of  Er;LiYF4  to 
undergo  self  Q-switching  on  the  transition  at  551  nm  [8]  when  pumped  by  multi-photon  absorption 
upconversicMi  processes  was  avoided.  Average  output  power  was  low,  only  2  mW  compared  to  20 
mW  obtained  in  cw  operation  under  identical  conditions  without  the  piezo-acoustic  modulator. 
Evidently  this  was  due  to  unexpectedly  high  losses  in  the  mode-locker.  Work  in  progress  on  a  new 
MgO:LiNbC)3  modulator,  designed  to  avoid  photorefractive  losses  caused  by  intense,  visible  intra- 
cavity  radiation  in  this  upconversion  laser,  should  lead  to  substantially  higher  output  power. 

In  the  second  part  of  this  work,  examination  of  the  upconversion  excited  by  1.06  pm  irradiation  of 
Tm-doped  fibCTS  revealed  the  dependence  on  input  power  shown  in  Figure  3.  Emission  intensity  in 
this  figure  is  directly  proportion^  to  1.06  pm  absorption  from  the  state  of  Tm^and  clearly 
indicates  that  absorption  increases  non-linearly  at  this  wavelength  as  excitation  is  increased.  The  figure 
also  gives  preliminary  evidence  of  a  threshold  for  this  behavior  at  very  low  intensity. 


Figure  2.  Mode-locked  pulse  width  versus  position  of  the  output  coupler.  The  solid  curve  is  the  best 
fit  of  a  quadratic  regression. 


0  Figure  t  Fluorescence  intensity  at  480  nm  ( >  G4  -» ^He)  from  Tm3+  impurities  excited  at  Xex=  1 .06 

|im  with  a  NdrYAG  laser  focus^  into  a  silica  fiber.  Core  diameter  of  the  fiber  was  5.1  |xm  and  dopant 
concentration  was  20,000  ppm.  The  core/cladding  index  difference  was  0.00675. 

DISCUSSION 

•  The  results  of  Figtfres  1  and  2  were  obtained  by  actively  mode-locking  a  gain  medium  with  an 
inversion  sustained  entirely  by  upconversion  of  infrared  radiation.  They  demonstrate  the  viability  of 
this  novel  approach  to  direct  generation  of  short  pulses  in  the  visible  spectral  range  from  solid  state 
sources.  Significant  improvements  can  be  anticipated  from  the  use  of  modulators  which  avoid 
photorefractive  damage  and  from  the  identification  of  room  temperature  upconversion  materials. 

#  In  addition  to  providing  stimulated  emission  at  new  frequencies,  cooperative  upconversion  processes 
also  give  rise  to  optical  nonlinearities  which  may  be  useful  in  novel  schemes  for  passive  mode-locking. 


For  instance,  the  data  in  Fig.  3  reveal  nonlinear  absorption  with  an  apparent  threshold  of  20  mW  and 
quadratic  intensity  dependence  on  incident  power,  indicating  that  such  nonlinearities  may  even  be 
useful  in  heavily-doped  fiber  lasers.  Saturation  of  an  ordinary  3-photon  ESA  process  seems  an 
unlikely  mechanism  for  Fig. 3  data  at  the  power  levels  of  this  experiment.  We  tentatively  conclude  that 
avalanche  absorption,  a  cooperative  effect,  accounts  for  this  interesting  behavior  at  the  Nd^-  laser 
wavelength.  While  three  photons  are  required  to  reach  the  'G4  state  in  principle,  upconversion  should 
indeed  be  quadratic  in  the  presence  of  an  avalanche,  since  the  onset  of  strong  cross  relaxation 
effectively  saturates  the  first  optical  transition  [6].  Drawing  on  these  results,  and  on  the  recent  repon 
of  insubilities  mediated  by  delocalization  of  excitation  in  upconversion  lasers  [9],  we  now  discuss  new 
ways  in  which  cooperative  nonlinearities  in  laser  gain  media  may  furnish  passive  mode-locking. 

First,  with  their  unusual  propeny  of  increased  absorption  with  increased  pumping,  avalanche 
processes  can  be  expected  to  suppress  cw  laser  action  and  promote  short  pulse  generation  in  a  simple 
way.  One  might  consider  for  instance  a  Tm,Nd  co-doped  fiber  laser  in  which  Tm  impurities  serve  to 
introduce  significant  losses  (above  avalanche  threshold)  for  cw  laser  operation  at  1.06  |im,  restricting 
operation  to  a  short  pulse  regime  characterized  by  saturation  of  the  ESA  at  1.06  |im.  In  such  an 
application,  the  repetition  rate  need  not  be  limited  by  the  avalanche  decay  time  but  rather  by  the  ESA 
saturation  recovery  time  which  is  typically  much  shorter. 

Internal  gain  suppression  at  1.06  pm  by  cooperative  absorption  of  Tm  could  also  be  used  to  promote- 
laser  operation  (cw  or  mode-locked)  at  alternative  Nd  wavelengths  with  lower  gain  cross  sections, 
such  as  0.9  or  1.3  pm.  At  these  wavelengths,  ellipse  rotation  caused  by  (off-resonant)  avalanche- 
induced  birefringence  would  be  experienced  by  intra-cavity  light  [10].  This  effect  could  then  mediate 
self  mode-locking  in  "nonlinear  mirror"  fiber  laser  designs  [11],  avoiding  avalanche  absorptive  losses 
at  the  lasing  wavelength. 

Finally,  pulsing  instabilities  due  to  delocalization  of  rare  earth  excitation  in  cooperative  upconversion 
Ipera  have  recently  been  reported  [9].  Their  repetition  rate  is  dependent  on  the  upconversion  rate, 
lifetimes  of  the  active  ion,  and  the  photon  density.  Population  pulsations  due  to  this  mechanism  have 
been  predicted  to  occur  even  in  upconversion  media  without  feedback.  Hence  instabilities  of  this  type 
are  self-starting  in  principle  and  have  a  calculable  prima^  frequency.  A  self-starting,  cooperatively 
mode-locked  laser  could  take  advantage  of  this  mechanism  to  provide  short  wavelength  sources.  Such 
a  device  might  be  realized  by  matching  the  cavity  round  trip  time  to  the  pulsation  period. 
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